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Abstract 

We describe the use of generalized unitarity for the construction of correlation func- 
tions of local gauge-invariant operators in general quantum field theories and illustrate 
this method with several calculations in M = 4 super- Yang-Mills theory involving BPS 
and non-BPS operators. Form factors of gauge-invariant operators and their multi- 
operator generalization play an important role in our construction. We discuss various 
symmetries of the momentum space presentation of correlation functions, which is nat- 
ural in this framework and give examples involving non-BPS and any number of BPS 
operators. We also discuss the calculation of correlators describing the energy flow in 
scattering processes as well as the construction of the effective action of a background 
gravitational field. 
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1 Introduction and summary 



Correlation functions of gauge-invariant operators are natural observables in both conformal 
and non-conformal field theories. In the early days of the AdS/CFT correspondence corre- 
lation functions of BPS operators played an instrumental role in establishing and testing it. 
More recently, correlation functions have been shown to exhibit fascinating relations to other 
quantities: limits in which the gauge invariant operators are null separated are related to 
expectation values of null polygonal Wilson loops with or without local operators [H [2j [31 E] 
and, in M = 4 super- Yang-Mills theory in the planar limit, the same limit is related to scat- 
tering amplitudes [5]. Moreover, since any curve may be approximated to arbitrary precision 
by a null polygon, it is in principle possible (though perhaps difficult in practice) that special 
limits of correlation functions with arbitrarily many operators can be used to construct the 
expectation value of Wilson loops of arbitrary shape. 

Correlation functions of gauge-invariant operators also have a natural place in non- 
conformal theories, such as QCD, where e.g. suitable correlation functions capture certain 
inclusive properties of the final state of scattering processes. Among them are the energy 
correlators, originally introduced in [B]. As discussed in [TJ, they have a counterpart in 
conformal field theories, where they describe the final state produced by the time evolution 
of some localized excitation and may be thought of as particular analytic continuations of 
regular correlation functions. 

Even though in conformal field theories their axioms guarantee that higher-point corre- 
lation functions are determined by the two- and three-point functions, an explicit evaluation 
along this line is not straightforward. It is therefore interesting to devise methods to directly 
evaluate them - at weak and at strong coupling - both for generic positions of operators or 
directly in special limits. 

The classic approach to the calculation of correlation functions makes use of Feynman 
diagrams in momentum space, position space or superspace; in M = 4 sYM theory three- 
point functions of scalar operators of dimension A < 5 have been evaluated systematically 
using such methods in [5J. An efficient reorganization of this approach is the Lagrangian 
insertion formalism [91 [10j [11]. In this theory, and in general in all theories in which the 
planar dilatation operator defines an integrable Hamiltonian, the calculation of three-point 
(as well as that of higher-point) functions benefits from use of integrable model techniques 
[T2"l [T31 dU [112 dS] . More recently it was proposed [TTJ and illustrated for the correlation 
function of four BPS operators in M = 4 sYM theory that light-cone superspace can be 
efficiently used for this purpose. The correlation functions of chiral stress tensor multiplets 
in M = 4 sYM theory enjoy special properties at the integrand level [18] - a permutation 
symmetry which becomes manifest when the integrand is constructed in the Lagrangian 
insertion formalism. This was used to devise an efficient method for the construction of 
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four-point correlation functions of stress tensor multiplets [TH] which led to the evaluation 
of the five-loop correction to the anomalous dimension of the Konishi multiplet ^2U\. The 
result confirms the integrability-based predictions [2T| [22]. 

In conformal field theories that have a string theory dual, certain classes of correlation 
functions of local operators may be evaluated at strong coupling using semiclassical expansion 
on the string theory side [231 121 ESJ [261 [23 EH EHl EQ] . It was moreover argued [HI [32] that 
AdS supergravity scattering amplitudes with suitable boundary conditions follow an on-shell 
recursion relation of the same type as the gauge theory on-shell recursion relations. 

Independently, increasingly efficient perturbative computational techniques - generalized 
unitarity, on-shell recursion relations, etc. - led to remarkable progress in our understanding 
of scattering amplitudes of M = 4 sYM theory and to the discovery of new and powerful 
symmetries - dual super conformal symmetry [3"3"l I3~4] . color/kinematics duality [35J - which 
severely constrain the scattering matrix. As we will explain in the next section, correlation 
functions of local gauge invariant operators may be interpreted as special scattering am- 
plitudes of sources in the theory obtained by adding to the action the operators and their 
sources. One might naturally expect that the modern techniques developed for S-matrix 
calculations can also be used to calculate correlation functions; as we shall see, this is in- 
deed the case. While correlation functions are naturally functions of the positions of the 
operators, an approach that mirrors the calculation of scattering amplitudes will yield their 
Fourier-transformed expressions, i.e. the momentum space correlation functions. Since po- 
sition space correlation functions are conformally invariant the momentum space correlators 
should also have this property (albeit non-manifestly) and thus should be annihilated by the 
Fourier-transform of the conformal group generators. 

The relation between correlation functions and scattering amplitudes implies that the 
conformal symmetry of the former becomes - in the null separation limit - the dual conformal 
symmetry of the latter [5] . Momentum space expressions for correlation functions may also be 
used to search for additional symmetries, which are hidden in the position space expressions^] 

Apart from constructing position-space correlation functions, momentum space correla- 
tors can also be used to construct quantum effective actions in background fields. As we 
shall briefly discuss in § [51 coupling an action with external fields is essentially equivalent 
to deforming the action by various operators, with the background fields acting as sources. 
Compared to the deformations relevant for the calculation of correlation functions, the only 
difference is that, depending on the specifics of background fields, sources may appear non- 
linearly. Nevertheless, scattering amplitudes of background fields, perhaps expanded for 

3 The relation to amplitudes might suggest a relation between momentum space symmetries of correlation 
functions and position space symmetries of amplitudes. Such a relation may however be obscured by the 
null limit, which is not very transparent in momentum space. While such a relation can exist only in the 
null limit, it is not clear why momentum space correlation functions cannot have additional symmetries. 
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small moment a/ weakly varying external fields, yield the terms relevant for the construction 
of their off-shell effective action. 

Recently, the color/kinematics duality [22] has emerged as an important property of 
color-dressed amplitudes in certain gauge theories with only adjoint fields and antisymmetric 
structure constant couplings. It states that, given some amplitude, there exists a presentation 
of its integrand such that, if the color factors of some integrals obey the Jacobi identity, then 
the numerator factors of those integrals also obey a Jacobi-like identity. There is by now 
substantial evidence in favor of the duality, both at tree level [36], EH EH EH] HQ] HI] and at 
loop level [12] . A consequence of this duality is the existence of nontrivial relations between 
the color-ordered partial tree amplitudes of gauge theory [35J, which have been proven both 
from field theory [13] and string theory [H] perspectives. We will explore whether a duality 
of this type exists for correlation functions. A natural expectation is that, if present, it would 
be most easily visible in momentum space correlation functions. 

The structure of the operators makes it unlikely that color/kinematics duality always 
holds for all internal edges of the integrals that appear. Indeed, color /kinematics duality 
acts naturally on integrals associated to graphs constructed from 3-point vertices. Due to 
the operator insertions however, correlation functions naturally have multi-point vertices. 
One may attempt to define the duality by resolving the higher-point vertices into sequences 
of three-point vertices in color space. In general however, these vertices carry both antisym- 
metric structure constants as well as the symmetric ones, d a b c , and their generalizations. It 
is therefore not clear whether for a generic correlation function Jacobi relations can exist 
for all internal edges. A notable exception] which appears consistent with color/kinematics 
duality is provided by a four-gluon operator with the color structure given by d a bcd- If a 
Jacobi identity is not present and color and kinematic factors are no longer linked, it is still 
possible that numerator factors of various integrals are nonetheless nontrivially related to 
each otheijE Regardless of these details, we will argue that color/kinematics duality should 
hold for all internal edges not directly connected to one of the operators up to contact terms 
collapsing at least one of these latter edges. 

In this paper we discuss the use of generalized unitarity for the construction of momen- 
tum space correlation functions of BPS and non-BPS operators. While we will be mainly 
concerned with describing the details and subtleties of this approach, we also illustrate this 
method by recovering the known example of the four-point BPS operators and also by con- 
structing infinite classes of new ones: such as the n— point function of BPS operators and the 

4 We thank H. Johansson for providing this example, discussions on this point and for sharing his insight 
into [45] . 

5 This is reminiscent of the /3-deformed N = 4 sYM theory where tree- level kinematic factors are 
related despite absence of a relation for the corresponding color factors. The importance of such relations 
remains an open question. 
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n— point function of BPS operators and some number of twist-2 non-BPS operators at the 
next-to-leading order. In our examples we will focus on the M = 4 sYM theory; as for scat- 
tering amplitudes however, the method we describe here may be applied with suitable care 
to all quantum field theories. Essential ingredients in this construction are the (super-)form- 
factors of local gauge-invariant operators as well as their generalizations involving several 
operators. Unitarity and generalized unitarity have already been applied to the construction 
of higher-loop form factors in |47J HSJ HH1 \5U\ . 

In § [2] we will describe the similarities between scattering amplitudes and correlation 
functions and the use of generalized unitarity for the construction of the latter for general 
gauge-invariant operators, the relevance of generalized form factors and the need for regular- 
ization and renormalization. We will also identify the components of correlation functions 
which may exhibit color/kinematics duality as well as may contain some hidden consequences 
of dual conformal invariance. In §[3] we collect the known expression of the MHV super- form 
factor of the chiral stress tensor multiplet and list (while relegating the details to appen- 
dices) the generalized (two-chiral stress tensor multiplet) form factor as well as the MHV 
super-form factor of scalar non-BPS operators and of the general twist-2 non-BPS operators. 
We will use them in the examples we discuss in § HI 

We construct the leading order and the next-to-leading order correlation function of 
four BPS operators in momentum space, Fourier-transform it to position space and recover 
the known results (511 [52]. While the calculations are mainly carried out in components, 
in § 14.1.31 we illustrate the use of manifestly supersymmetric methods for the construction 
of correlation functions. We also derive the (connected part of the) next-to-leading order 
correlation function of n BPS operators and check that its null limit reproduces the n-point 
MHV amplitude, as originally shown in [5]. In § 14.21 we compute (the connected part of) 
the correlation function of one and two general non-BPS twist-2 operators and an arbitrary 
number of BPS operators. Similarly to the discussion of correlators of BPS operators, we 
verify explicitly that their null limits reproduce the n-point MHV amplitude. We also discuss 
the structure of the correlator of m twist-2 and n BPS operators in a split configuration. 
For the three-point function ^ we discuss the appearance of the anomalous dimension of the 
non-BPS operator, its renormalization, as well as the vanishing of the three-point function 
in limits in which the twist-2 operator becomes a conformal descendant of a BPS operator. 

In § [5] we include a general discussion of correlation functions involving stress tensors, 
which may be used to determine the effective action in non-dynamical gravitational back- 
ground. We also discuss the energy correlators, which capture the energy flow in the time 
evolution of some composite state. After recalling their definition, we illustrate this non- 
standard type of correlation functions by evaluating several examples to first nontrivial order: 

6 While writing up this paper we received 53 in which the three-point function of one twist-2 non-BPS 
operator and two BPS operators was computed through a different method. 
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the one- and two-point energy correlator in the state created by a BPS scalar operator as 
well as the one- and two-point energy correlator in the state created by a twist-2 non-BPS 
operator. Section § [6] contains our conclusions and a summary of our results. Several ap- 
pendices contain detailed derivation of some of the formulae used in the main body of the 
paper. 



2 Generalized unitarity and correlation functions 

Generalized unitarity has been developed as a very efficient tool for the calculation of on-shell 
scattering amplitudes in quantum field theories. From a modern perspective it is interpreted 
as a specific organization of the Feynman graphs contributing to an amplitude which exposes 
the vast simplifications occurring when off-shell Green's functions are amputated and placed 
on shellfl each generalized cut is, on the one hand, a product of tree-level amplitudes and 
on the other it is the subset of Feynman graphs that contain the cut propagators. From this 
perspective it appears natural that anything that is expressible in terms of Feynman graphs 
and is gauge invariant may be constructed through generalized unitarity-type methods^] 
An example in this direction are form factors; particular form factors have recently been 
considered at various loop orders in [3TJ HHJ EU ESI EDI EH] in Af = 4 sYM theory as well as 
in QCD |56l E7]. 

While this interpretation of generalized unitarity makes clear its applicability, it is not 
difficult to give form factors the interpretation of scattering amplitudes. In a similar spirit 
correlation functions, which are the focus of our paper, can be given the same interpretation.^] 



2.1 Correlation functions as scattering amplitudes 

Let us start with the path integral expression of an n-point correlation function in some 
quantum field theory defined by the Euclidian action 5s[$], 

(0 1 (x 1 ) . . . O n (x n )) = J [£>$] O^arx) . . . O n (x n )e~ s ^ , (1) 

where $ generically denote the fields of the theory and the local operators Oi are gauge- 
invariant combinations of them. All correlation functions may be packaged into a generating 

In the unitarity construction of scattering amplitudes based on the optical theorem the imaginary part 
of an amplitude has this property manifestly; the real part is constructed through a dispersion integral. 

8 It is presumably possible to generalize this statement to gauge- variant quantities at the expense of having 
non-vanishing contributions from ghosts and unphysical degrees of freedom. 

9 From the perspective of the optical theorem, this implies that correlation functions should also be 
expressible in terms of dispersion integrals. 
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functional, by introducing local sources for all possible/desired local operators: 

Z[J U J 2 ,...\ = J [£>$] e-^ffl-JAE^iW^W . (2 ) 

Correlation functions of specific operators are then extracted by differentiating with respect 
to the relevant sources and setting all sources to zero: 

(0 1 {x l )...O n {x n ))= 5 \ S Z[J U ...\ . (3) 

0jl{Xi) . . .bj n (x n ) Ji-+0 

The sources Ji may be interpreted as non-dynamical fields. If Se is the action of a conformal 
field theory with a closed string theory dual [SB] , then the sources may be interpreted as the 
boundary values of the fields describing the closed string states [591 EH] • In particular, for 
BPS operators, they are just the boundary values of the fields of the relevant supergravity 
theory. 

The construction above is formally identical to that of scattering amplitudes; the only 
difference is that, for scattering amplitudes, J are sources for the fundamental fieldsB From 
this perspective, if we interpret the sources J as fields, we may also interpret the correlation 
functions of operators as scattering amplitudes of their associated source-fields J in the 
theory with the modified action 

S%° d [$,J] = S E m + [d d xJ2 Ji{x)Oi(x) . (4) 

We may therefore view the scattering amplitudes of source-fields Ji as the correlation func- 
tions of the corresponding operators Oi. If the field theory has a string theory dual we may 
interpret this as the scattering amplitude of closed strings with specified boundary condi- 
tions. Similarly, we may interpret the scattering amplitudes of one source-field J and any 
number of fundamental fields $ as the form factor of the operator corresponding to that 
source, 

(0|$i---$ m ), (5) 

as also mentioned in [55]. From the perspective of a string theory dual we may interpret 
this as the scattering amplitude of a closed string into open strings. Amplitudes with several 
sources and fundamental fields have a similar interpretation in the context of gauge/string 
duality; from a field theory perspective we will refer to them as "generalized form factors" 
or as "multi-operator form factors". 



10 Also, scattering amplitudes are constructed in Lorentzian signature with the time-ordered operators. 
For correlators in Euclidian signature one may formally introduce a radial ordering of operators. 
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It is important to stress that these are not necessarily related to form factors of multi-trace operators. 
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The gauge invariance of the operators Oi and of the corresponding source fields Ji guar- 
antees that this action is BRST invariant. As in the absence of source fields, it is possible 
to use a sequence of Ward identities [61] to show that scattering amplitudes with external 
unphysical fields or ghosts do not contribute to unitarity cuts, 

As described above, while J are interpreted as fields, they are nevertheless non-propagating 
and thus they may appear only as external lines of an amplitude and do not have a well- 
defined on-shell condition. To make an even closer analogy with scattering amplitudes one 
may formally promote them to propagating fields by assigning them massive quadratic terms 
with different masses for each source; they may then appear both as internal and external 
lines. Restricting the integrands of the resulting amplitudes to terms that have no poles 
as the mass of the sources is varied guarantees that all contributions with internal sources 
are projected out. Last, the remaining mass dependence is solely associated to the norm of 
the external momenta and should participate in the Fourier-transform to position space. 



2.2 On the presentation, structure and symmetries of correlation functions 

Correlation functions of operators with definite dimension transform covariantly under po- 
sition space conformal transformations and are strongly constrained by them. However, 
as for scattering amplitudes, the generalized unitarity method yields scattering amplitudes 
for source fields carrying definite momenta and thus gives the momentum space form of 
correlation functions 

(OM...OM) (6) 

where Oi(qi) are the Fourier-transform of the usual position space operators 

Position-space conformal invariance will be hidden in momentum space; nevertheless, since 
momentum space correlation functions should not have holomorphic anomalies^] [63] (unlike 

12 Thc required Ward identities do not hold if at least one of the operators - and hence the deformed action 
(|4j) — is gauge- variant. 

13 An example in this direction is the interaction term HTr[F 2 ] generated at one- loop level through a 
top-quark loop in the Standard Model (and it is responsible for the gluon-fusion production of the Higgs 
boson). The two- loop amplitude H — > ggg in this theory, i.e. the two- loop form factor of Tr[F 2 ] with three 
external states, was discussed recently in [SSJ [SU [57]- By formally taking the Higgs mass to infinity one 
forces the Higgs boson to appear only as an external state. The resulting amplitudes are either generalized 
(in the sense described below) form factors of Tr[F 2 ] or correlation functions of these operators. 

14 Holomorphic anomalies are likely to appear when the operator's momenta are restricted to be null. 
Since position space correlation functions should be conformally invariant, we expect that the constrained 
momenta should not play an important role in the Fourier-transform. 
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scattering amplitudes of fundamental fields) one may test whether they are annihilated by 
the momentum space form of conformal generators - in particular the special conformal 
generator 

n 

J2 K i ( 8 ) 

i=l 

d ( d 2 ■ d d \ d 

-2(A, - 3)— - qt^-K~ v ~ H^i^T + 2^(2 V 

where the sum runs over all operators, A is the conformal dimension and S l represents 
SL(2, C) on the %— th operator. Conformal symmetry should emerge as a manifest symmetry 
upon inverse Fourier-transform to position space, if the operators are chosen to have definite 
dimension. 

The unitarity method provides an efficient framework for systematically constructing 
and verifying the expression for any multi-loop amplitude in a massless field theory. This 
method, along with various refinements, has already been described in some detail elsewhere 
[64l l65l l66| 1671 l68l [69] ; see also [70j [TTJ [72] for recent reviews. Here we will discuss the 
additional information needed for the construction of correlation functions of local operators. 
A color-dressed generalized unitarity cut is a sum over products of color-dressed amplitudes 

iC Ai) A m---Am) (9) 

states 

where the c cut lines are placed on shell and we included a factor of % for each cut propagator; 
each cut line appears twice, leaving one amplitude factor and entering another. While not 
necessary, the amplitude factors Au\ can be chosen to be at tree level. In our case of cuts of 
scattering amplitudes of source fields, the factors Au) are either amplitudes of fundamental 
fields or (generalized) form factors. The cut construction of an amplitude formally proceeds 
by matching all cuts onto an ansatz in terms of Feynman integrals. 

A systematic strategy, which is designed to keep under control the size of the ansatz, is 
the maximal cut method in which one begins by first constraining the ansatz to reproduce 
all the (3L + (n — 3))-particle cuts (i.e. the maximal cuts, in which the maximal number of 
propagators for an n-point amplitude are cut) and then systematically proceeds to relax the 
cut condition on one propagator (next-to- maximal cuts), two propagators (next-to- next-to- 
maximal cuts) and so on. At each step one may reduce the expression of the cuts to cuts of 
master integrals by generalizing the methods of refs. [73] 157] beyond two-loop order and to 
all massive external legs. The advantage of this method is that, for any one cut, only a small 
part of the ansatz is relevant. For the resulting expression to be correct it must reproduce 
all generalized cuts. Since some cuts are special cases of others, it suffices to verify that 
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it reproduces a spanning set of cuts — i.e. a set of cuts that guarantee that all the other 
ones are satisfied. For simple (low-order) cases one may construct amplitudes by analyzing 
directly the relevant spanning set of cuts. 

It is therefore clear that an essential ingredient in the construction of the (momentum 
space form of) correlation functions are the form factors and generalized form factors of 
the corresponding operators. While it is obvious that there exists a very close relation 
between off-shell Green's functions of fundamental fields and correlation functions of local 
operators, the generalized unitarity makes it clear that a similar though slightly weaker 
relation exists between scattering amplitudes of fundamental fields and correlation functions. 
Let us consider the generalized cut shown in fig. [TJ the external lines are attached to blobs 
representing color-dressed form factors while Feynman diagrammatics guarantees that the 
blob with no external lines is some multi-loop scattering amplitude of fundamental fields. 
This structure is independent of the number of external legs of each form factor, with the 
highest-loop amplitude appearing together with the form factors with the smallest number 
of external legs. Other cuts are necessary in order to identify the terms in which propagators 
exposed in fig. [T] are collapsed. It is moreover necessary to consider cuts involving generalized 
form factors. The terms captured only by these cuts are those in which there is no (sequence 
of) propagator(s) between two operators; upon Fourier-transform to position space, such 
terms lead to contact terms which we know are relevant only in the OPE limit; it is possible 
that at least some contact terms may be inferred from the requirement that the position 
space Euclidean correlator satisfies all relevant Ward identities (see e.g. [H] for a discussion 
of 3-point functions and the determination of contact terms and [73] for the contact terms 
required by conformal invariance). 




Figure 1: Generalized unitarity guarantees the presence of multi-loop amplitudes in the 
correlation function calculation. Additional contributions involves at least one form factor 
with more fields than the corresponding operator. 
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This structure of correlation functions points to the fact that, in theories whose color- 
dressed scattering amplitudes of fundamental fields exhibit larger symmetries, such as M = 4 
sYM theory, momentum space correlation functions should be at least partly constrained by 
these symmetries. Among the remarkable properties of M = 4 sYM theory are the dual 
conformal invariance and color/kinematics duality of its scattering amplitudes. They may 
have some consequences on correlation functions as well. Indeed, the cut in fig. [T] in which 
the form factors have the smallest number of external lines implies that there should exist 
a presentation of correlation functions such that the numerator factors of the contributing 
Feynman integrals obey Jacobi identities for all internal lines not attached to an operator, up 
to contact terms that collapse one of these lines. Depending on the structure of the operators, 
in the complete correlation function the internal legs directly attached to any one of them 
may not necessarily obey a Jacobi identity and a case by case study appears necessary. For 
example, it was shown [76] that scattering amplitudes with an insertion of Tr[F 3 ] operator 
- i.e. the zero-momentum form factors of this operator - exhibit color/kinematics duality; 
the arguments above imply that we should expect that the correlation functions of several 
such operators should also exhibit color /kinematics duality in the zero-momentum limit. 
The fate of color/kinematics duality for the correlation function of Tr[F 3 ] operators with 
generic momenta is an interesting question, albeit one which we will not explore here. In 
general, color Jacobi identities exist only in the special case in which the color factors of 
operators contain at most one factor of the symmetric structure constants or their higher-rank 
generalizations and thus the correlation functions of only such operators may be expected 
to obey color/kinematics duality for all internal lines; the possible appearance of contact 
terms in Jacobi transformations on internal lines requires nevertheless a case by case study 
to ascertain whether color/kinematics duality is in fact realized. This structure is relatively 
similar to the properties of amplitudes in the (3— deformed M = 4 sYM theory [46] . 

Similarly, if excising the operators leads to a planar amplitudd^l. then the terms in the 
correlation function which contribute to this cut may exhibit some trace of the dual conformal 
invariance of that amplitude. In the complete correlator the non-planar terms (if present) as 
well as the presence of the operator insertions will break the usual dual conformal invariance. 

Regularization and renormalization are important issues which need to be carefully con- 
sidered. Since all external lines of momentum space correlation functions (or, equivalently, 
of source field amplitudes) are effectively massive, no infrared divergences can appear and 
thus infrared regularization is not necessary. Ultraviolet divergences have two possible ori- 
gins: divergences due to the structure of the undeformed Lagrangian and divergences due 

15 It is important to note that, since the source fields J are color-singlets, planar integrals are only a 
subset of the leading terms of the larger N c expansion. Indeed, if the non-planarity of a integral arises 
solely because external legs (or sources) are attached to internal lines, then that diagram yields, in fact, 
leading-color contributions. This has already appeared in the calculation of higher- loop form factors [481150] . 
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to the presence of the deformation by the gauge-invariant operators. The former may be 
eliminated by standard renormalization of the Lagrangian and lead to /3-functions for the 
various coupling constants. In M = 4 sYM theory they are absent. The later are eliminated 
by renormalizing the deformations and are solely related to the anomalous dimensions of 
operators. In particular, if one used renormalized operators, with the renormalization fac- 
tors given in dimensional regularization by the usual expression in terms of the anomalous 
dimension of the operator 



(g is the gauge coupling constant), the correlation functions would be finite in the ultraviolet 
as well. In general the anomalous dimension 7 requires a separate calculation. In M = 4 
sYM in the planar limit they are provided by integrability. 

One therefore has several options for accounting for the required renormalization: (1) 
one uses from the outset dimensional regularization and renormalized operators or (2) one 
carries out the calculation in four dimensions, regularizes dimensionally (if necessary) the 
integrals appearing in the final result and searches for potentially missing /i-type integral 
such that multiplication by the appropriate Z-factors (jTUl) renders it finite as e — > or (3) 
one uses another regulator, such as a higher-derivative regulator, for both the calculation of 
the amplitude and the calculation of the renormalization factors. 

While the strategy outlined here applies to generic quantum field theories, in the follow- 
ing we will restrict ourselves to M = 4 sYM theory. In this case the sources break maximal 
supersymmetry to the subalgebra that leaves invariant the deformation (T4J). BPS operators 
preserve some amount of supersymmetry, their anomalous dimension vanishes identically 
and their Z-factors equal unity. In their case no regularization is necessary. For the cal- 
culation of correlation functions that include non-BPS operators the deformed action (J3J) is 
formally nonsupersymmetrij^l. Nevertheless, since supersymmetry breaking is confined to 
the operator insertions, we expect that most of the features of maximally supersymmetric 
calculations will continue to exist here as well. Since to leading order (tree-level in position 
space) all correlation functions are rational functions and at the next-to-leading order only 
one-loop bubble integrals are divergent, a systematic regularization procedure starts being 
necessary only at the next-to-next-to-leading. 

Quite generally, operators with definite anomalous dimension are (complicated) linear 
combinations of simpler single-term operators with coefficients given |7?] , in the planar limit 

16 We recall here that ^t-type integrals are integrals whose integrand vanishes identically when evaluated 
in four dimensions. 

17 The same is true, in fact, for the action deformed by several BPS operators each of which preserves a 
different non-overlapping or partly overlapping subsets of supercharges. 




(10) 
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and for operators with at least one large charge, in terms of the solution to the Bethe 
equations [7S]. The form factors of the former operators are linear combinations of the form 
factors of the latter operators with the same coefficient^. To compute correlation functions 
of such operators we evaluate the correlation functions of the generic terms in each of them 
and then sum them with the appropriate integrability-determined coefficients. This strategy 
applies both to 3-point as well as to higher-point functions. 

As noted in [SI], scattering amplitudes may be interpreted as the (generalized) form 
factors of the zero-momentum Lagrangian. Similarly, form factors with more external legs 
than fields in the operator may be interpreted (e.g. through the MHV vertex expansion) as 
generalized form factors of the operator and a suitable number of zero-momentum on-shell 
Lagrangians. In particular, scattering amplitudes with more than four external fields can be 
interpreted as generalized form factors of several on-shell Lagrangians. It then follows that 
the generalized unitarity-based construction of the integrand of the next fc -to-leading order 
momentum space correlation functions of some operators is the same as the construction 
of the leading order correlation function of those operators and k on-shell LagrangianJ^I. 
We therefore see a formal parallel between the generalized unitarity calculation and the 
Lagrangian insertion method of [HI EU E] • 

The various monomials that enter the expression of operators with definite anomalous 
dimension may not always have the same number of fields. A simple example is the stress 
tensor, whose expression contains a bilinear in the field strength which has two, three and 
four fields. This is, in fact, the generic structure of terms on the higher levels of supersym- 
metry multiplets with the increase in the number of fields being due to non-linear terms 
in supersymmetry transformations. Clearly, these terms contribute to form factors with 
different number of external legs; in the standard organization of the Lagrangian, in which 
the coupling; constant appears as an overall factor, these terms contribute to different loop 
orders, cj It has been suggested in [7H1 EDI E] their contribution cancels against contact 
terms arising from terms proportional to the equations of motion from Lagrangian insertions 
and, at least for the purpose of constructing the null limit of correlation functions, they 
may be set them aside. Their close relation to contact terms (which in a momentum space 
framework should appear as Feynman integrals with cancelled propagators) suggests that a 
shortcut to determining the contribution of nonlinear terms in supersymmetry transforma- 
tions to correlation functions may be imposing the Ward identities of the various symmetries 
on (the Euclidian) correlation functions found in the absence of these terms [7H [TSJ [11] . 

18 Since the L-loop mixing coefficients arc proportional to A i_1 , the L-loop form factor requires use of the 
(L + l)-loop eigenvectors of the dilatation operator. 

19 At each order beyond the leading order one needs to add one additional internal line and hence one 
additional operator with at least three fields. 

From a position space perspective one might be tempted to interpret all of them as contributing to the 
same order - e.g. at tree level - as they do not involve any integration. 
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3 Some single-operator and multi-operator form factors 

We have argued in the previous section that tree-level form factors and generalized form 
factors are essential for the construction of correlation functions. They may be constructed 
in several ways, mirroring the various methods for the construction of tree-level scattering 
amplitudes. As in that case it is useful to assemble them into super-form factors; they are 
labeled by the coordinates of two superspaces: the usual on-shell M = 4 superspace, with 
Grassmann coordinates r\ A and an off-shell superspace for the multiplet of operators, with 
Grassmann coordinates 7. A convenient one is the harmonic superspace [82]; the on-shell 
superspace fieldg^] may also be written in harmonic superspace form [SOj [81]. Whenever 
necessary, we will follow the notation there as well as in [54"[ [S3] and collect some details in 
Appendix [A] 

A possible approach to the construction of (super-)form factors is to use the MHV vertex 
rules. This strategy was first applied to the construction of form factors of a particular gluon 
operator in [SI] and more recently in [85]. This method requires independent information on 
MHV form factors; they are defined to be those that have the minimal number of fermionic 
coordinates rj. This definition mirrors that of MHV amplitudes. 

BCFW-like recursion relations may also be used if a suitable shift can be found. In this 
approach only the tree-level form factor with the minimal number of fundamental fields is 
necessary. As in the case of scattering amplitudes of fundamental fields, different BCFW 
shifts yield different presentations of the same form factors. 

In this section we collect several examples of form factors and generalized form factors 
with arbitrary number of external states which will be useful in the examples of correlation 
functions we will discuss in the next section. 



3.1 The tree-level MHV form factors and generalized form factors of the chiral 
stress tensor multiplet 

The tree-level MHV super-form factor of the chiral part of the stress tensor multiplet 

T was found in ref. [5U E3] through a combination of symmetry constraints and explicit 
calculations: 

n 

i=l \ ' ' % i % 

21 Specific on-shell external fields may be extracted in the usual way, by specifying their position in the 
on-shell Af — 4 multiplet 

Hv) = 9- + VA^ A + l-.ri A riB<j} AB + ^VaVbVcc ABCD iPd + ^VAVBVcVDe ABCD 9+ ■ 
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The two 5— functions may be combined into a single one by multiplying their arguments with 
the suitable harmonic variables: 



i=l 



The form factor of the chiral primary operator is extracted as the coefficient of (7+) 4 while 
that of the on-shell Lagrangian as the coefficient of (7+) . 

The simplicity of the MHV form factor resembles that of MHV amplitudes. As discussed 
in [SI], the conjugate super- form factor may be obtained by conjugation and Grassmann 
Fourier transform of the super-form factor of the chiral primary operator: 

i 1 t=l L J 

/ \[d%^5^{Y,~ X ^t<A) ( 13 ) 

J i i 

An alternative expression, in terms of a single auxiliary Grassmann integral, is 



N max MHV 



n 1 r 

8\q - £ - E W II T-TTT / d '< 5(2) (^Af - ^) • (14) 

i i i=i L l > 1 + A J ./ 



One might wonder whether there exist generalized MHV form factors — i.e. MHV form 
factors with several operator insertions. The supersymmetry Ward identities discussed in 
[5^] may be easily extended to such cases: for the chiral stress tensor multiplet, a possible 
solution is obtained by simply replacing 7+ by the sum of the 7+ coordinates of all the 
operators, Y^il+'X- 

^» v r fc = F*%-E A ^)^ , (15) 

i j=l i i 



with some coefficient Fl 22 l. This expression however contains at most four 7+ coordinates and 
therefore will not contain e.g. the generalized form factor of two chiral primary operators 
(CPO-s) (which would require eight 7+ coordinates^ 

We may alternatively consider a generalization of eq. f lT2|) in which we replace l^ a 7j" by 
Ylk=i ^Aal%+' However, extracting various components - such as the two-CPO component - 
out of this expression one quickly finds that they cannot be generated by Feynman diagrams; 
for example, the two-CPO form factor will not have any scalar external lines. 



22 The supersymmetry Ward identities discussed in [54] suggest that F = rim=i( TO ' TO + -0 • 

23 It contains however the generalized form factor of the CPO and any number of chiral Lagrangians. 
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Last, one could potentially imagine an MHV form factor with different harmonic vari- 
ables for each operator that would not allow for a representation in terms of a single eight- 
dimensional Grassmann 5— function as it would only satisfy the supersymmetric Ward iden- 
tities up to equations of motion. However, simple counting suggests that the generalized 
form factor of two CPO-s cannot appear in an MHV generalized form factor because the 
relevant form factor should have at least twelve Grassmann S— functions (eight to saturate 
the eight 7 integrals isolating the two chiral primaries and four isolating two external scalar 
fields). 



3.2 Generalized tree-level NMHV form factor of two BPS operators 

Generalized form factors, as well as the N fc MHV form factors of (fc+1) BPS operators, are an 
important ingredient in the construction of higher-loop corrections to correlation functions 
of CPO-s; in particular, they are necessary to identify potentially degenerate contributions 
in which propagators between operators are canceled. The evaluation of this form factor 
turns out to be most efficient through the MHV vertex expansion [551 El] • We should sum 
over all graphs with two vertices in which either vertex is a super- form factor, fig. [2£ a), or 
one vertex is a generalized two-operator MHV form factor and the other one is a regular 
MHV amplitude, fig. EJ^b): 

-r-NMHV _ -r-NMHV , tNMHV 

j- — + y-( 6 ) . i^idj 

Since, as we argued in the previous section, an MHV generalized form factor for the lowest 
component of the chiral stress tensor multiplet does not exist, the diagrams of the second 
type have vanishing contributions 

j:NMHV = (17) 

For higher components of the chiral stress tensor multiplet this contribution may be nonva- 
nishing; we will however not be interested in it here. 

Relegating the details to Appendix El we list here the final result for J 7 ^^ 4 ^'- 

n 

jrNMHV = 5 8 (l^l + 2t 1 l + -J2v m X m )S\q i + q - 2 -p 1 ----p n ) 

m=l 

2 n i-2 2 n 

w=l i=l j=i+l w=l i=l 



3! lX=i( mm + 1) 



1(3 




Figure 2: The two generic terms in the MHV vertex expansion of the generalized super- 
form factor with two insertions of the chiral stress tensor. Grey blobs represent regular 
single-operator MHV super-form factors, the hashed blob represents an MHV generalized 
super-form factor and the white blob represents a regular MHV scattering super-amplitude. 
The momenta carried by the two operators are qi and q^- 



where Afj and Bf are defined as: 

4 = ^(vTimT - «7£+lOSl<] + E^(r|gg|i])K 6 ( 7 ^ + z) -j^VAr(n)) 

r=i r=i 

x «<^> - 1 v Br {n)) («&<•&.■> - ± ^) ^not^mm\) ■ (19) 

Q 7/ ,+ „.,+ b'c' TABC / 

Bf = M ' 2 ^ )2 ( - iZn^TiHtXT^i) - l e a ' d '<^^(7^7^)(7^7^ + ) 

+VTiw£ b e b c (^ + i)5 ad (^ + ^ + ) -VciW$ b e b c (^ + i)5 ad (^ + ^ + )^J . (20) 

It is not difficult to check that the simplest such form factor has two external scalars and it 
has no momentum dependence, as one might expect based on Feynman diagram calculation. 



3.3 Tree-level MHV form factor of a non-BPS scalar operator 

As discussed in § [2], one of the building blocks of correlation functions of non-BPS operators 
are their tree-level form factors and generalized form factors. Non-BPS operators with 
definite scaling dimension are typically linear combinations of simpler single-term operators 
with the same quantum numbers. Their relative coefficients are functions of the coupling 
constant and, at least for long operators, may be determined using the integrability of the 
dilatation operator of M = 4 sYM theory [771 EH] - The form factors of non-BPS operators 
may be therefore interpreted - order by order in perturbation theory - as the sum of form 
factors of these single-term operators with the same coupling-dependent coefficients. We 
thus need to focus on them. 
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The operator mixing determining the eigenvectors of the dilatation operator is, in general, 
rather complicated with the only constraints arising from charge conservation and the fact 
that such mixing can occur only between operators with the same classical dimension. There 
exist, however, "closed sectors" [87] in which operator mixing involves a very restricted class 
of operators. An example is the so-called S77(2)-sector which contains operators of any 
dimension constructed out of two complex scalar fields which are not conjugate of each 
other. 

The R-symmetry properties of the M = 4 sYM fields implies that, at tree- level, it is 
more efficient to evaluate the form factor of scalar operators of arbitrary charges rather than 
restricting ourselves to specific charges . We carry out this calculation in Appendix [B] using 
a BCFW recursion relation. Defining 




H-aAbB = VAaVBb ~ VBaVAb + ^abVBaVAb, (21) 

(EaibiTj = KiK^M ' ( 22 ) 

the MHV super-form factor of the operator OA 1 B 1 ...A k B k = Tr(0 j4lBl • • -4>A k B k ) with n exter- 
nal fields is given by: 

where the first sum runs over all over the sets {ai,6i--- ,cik,bk} with a\ < b\ < a2 < 
&2 • • ■ bk-i < a>k < bk and its cyclic permutations and Sp stands for the trace over the 
SL(2, C) indices of the product of E matrices. By restricting the pairs (Ai, Bi) to only two 
values, e.g. {(1,2), (1,3)} or {(1,2), (1,4)}, one finds the SU{2) sector operators and their 
form factors. 

The tree-level form factors ( 12"B"|) can be used to construct correlation functional; one 
first constructs the correlation functions of the single-term operators O A 1 B 1 ...A k B k - l then one 
takes appropriate linear combinations with integrability-determined coefficients and finally 
one renormalizes the result by multiplying with the corresponding factors Z 7i , cf. eq. ( fit)]) . 

3.4 Tree-level MHV form factor of a general twist-2 operator 

In the next section we will construct examples of n-point functions with one twist-2 spin-5* 
operator. These operators are linear combinations of O^s'x = Tr(.D^0 AB ''D^T 1 (j) CD ) 

Otr° = E <*,n{\)Oi$S D , (24) 

n=0 



24 Thcy could, of course, also be used to construct higher-loop form factors for these operators, determine 
their mixing and anomalous dimensions, etc. 
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where the coefficients c n (A) are determined by requiring that these operators have definite 
anomalous dimensions. At one loop and if (f> = <p CD = Z, they are given [88] in terms of 
the Gegenbauer polynomial^! 

2 , s = (n ■ D, + n ■ D 2 ) s C'J 2 (*g^W[Z(&)Z(fc)] , (25) 

v ' €1=6 

where n ■ D are covariant derivatives in the adjoint representation in the light-like direction 
specified by the vector n; one possible choice is n ■ D = D + . The one-loop coefficients in 
are cs !n ,o = ( — l) n ( n) 2 ; the two-loop coefficients may be found in 



We will describe in Appendix[B]the evaluation of the form factors of the operators O^g' 00 
through a BCFW recursion relation. It is convenient to express the result in terms of the 
^-dependent combinations H in eq. ff2T]) : 



■■■»> = ff'faE+i) E H — H — EE I E iv- 1 lEs 



c-1 a-1 / k \ x , , x S-x 



ni=iK m + i ) 



{a,b,c,d} k=b l=d \r=l+l / \s=k+l 



x ( (bW-p_k\c) + (b\p k+ i*-\c) _ {bc) \ f (d\a- pi \a) + (d\p l+1 a~\a) _ (&) } (2g) 

where q is the momentum of O, the sum runs over all sets {a, b, c, c?} where a < b < c < d 
or d < a < b < c etc. For small number of external particles this expression may be easily 
verified using Feynman graphs. As in the case of the scalar non-BPS operators, taking the 
appropriate linear coupling constant-dependent combinations of these expressions yields the 
form factors of twist-2 spin-S* operators with definite anomalous dimensions. 



4 Examples of correlation function construction 

Using the form factors constructed in the previous section we shall construct examples of 
correlation functions to leading order (LO) and next-to- leading order (NLO). Some of these 
correlators have been known for some time and we reproduce their expressions. We will also 
discuss various limits and properties of our results. 



4.1 Correlators of BPS operators 

Correlation functions of four chiral stress tensor multiplets have been evaluated to high loop 
order [191 120] using a hidden symmetry of their position space integrand [18] . We will discuss 

25 The Gegenbauer polynomials appearing for operators discussed here are the same as the Legendre 
polynomials with index S. 
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Figure 3: Disconnected and connected quadruple cuts. 



here, to a much lower loop order, a unitarity-based approach to the same correlator. As we 
will see, this construction generalizes quite easily to correlators of any number of operators. 

The simple structure of the CPO-s implies that the construction of their 4-point function 
through generalized unitarity is almost identical to that in terms of Feynman diagrams. 
From (fTTj) it follows that the (MHV or anti-MHV) form factor with two external lines of 
the chiral primary operator has no momentum dependence, which is the same as the vertex 
containing the source of the operator and two scalars. For more operators we need to use 
the generalized form factors; nevertheless, for scalar operators, they are quire similar to the 
results of Feynman diagram calculations. 

R-charge conservation implies that the correlation function of four CPO-s (two-index 
symmetric traceless dimension-2 operators) is determined by six functions of the positions 
of the operators [5Tj : 



(Tr[0 J V ]Tr[0 J2 j2 ]Tr[0 j3 j3 ]Tr[0 Y4 j4 ]) 

= 01^12^34 + a 2<$13^24 + «3^14^23 + ^1^13^14^23^24 + ^2^12^14^32^34 + &3 ^12^13 &42&13 



(27) 



with 



'{h\ u {J± u h} u \Jz} ' ( 28 ) 

where ij and Jj are SO (6) indices. As mentioned before, since the sources are singlets under 
the M = 4 gauge group, each term in the momentum space correlator may be extracted from 
un-ordered source-field scattering amplitudes. 



4.1.1 Leading order 



As it is well-known, at this order the contributions to the (Fourier-transform of the) coef- 
ficient functions a» come only from disconnected graphs. In our language, the generalized 
(quadruple) cut contributing to e.g. a\ is shown in fig. 3(a) The b{ coefficients receive con- 
tributions only from connected graphs. To determine them one may consider two-particle 
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cuts with the 2-operators form factors discussed in § 13. 2[ A close inspection of these form 
factors with two external fields implies however that there is always a propagator between 
the two operators; this in turn implies that the leading order correlator is determined by its 



quadruple cuts, such as the one shown in fig. 3(b) and its non-cyclic permutations. Depend 



ing on the specific choice of R-charge carried by the four operators only some of these cuts 
may exist. 

Choosing as representatives of the four CPO-s the operators 0\ = Tr[0 12 12 ], 2 = 
Tr[0 34 34 ], 3 = Tr[(p 13 (p 13 ] and C 4 = Tr[0 24 24 ] only the disconnected quadruple cut is 
non-zero. This correlation function determines the coefficient a±. The relevant form factors 
which determine the quadruple cut are extracted from ffTTj) (by appropriately choosing the 
harmonic variables) and are constant; it is therefore trivial to find that the momentum space 
representation of this correlator and thus of 03 is @ 

S 3 = (d 1 (g 1 )d 2 (g 2 )d s (g 3 )d i (g i )) (29) 

= S ^ + ^ + ^)J^^^I^ W ^- (30) 

It is easy to Fourier-transform this expression to position space before the p and I integrals 
are carried out, with the expected result. By permuting the labels of the various operators 
it is not difficult to find the coefficients a x and a 2 . 

Choosing as representatives of the four CPO-s the operators 0\ = Tr[<j) 12 <j) 23 ], 2 = 
Tr[0 14 24 ], O3 = Tr[(/> 13 </> 13 ] and O4 = Tr[0 24 </> 34 ] the only non-zero quadruple cut is the 



one shown in fig. 3(b) From the R-charge assignment it is easy to see that this correlator 
determines the coefficient b 2 in eq. (1271) . The quadruple cut may be easily evaluated to be 
consistent with 

4 

~b 2 = (d 1 (q 1 )d 2 (q 2 )d 3 (q 3 )d,(q 4 )) = 5 4 (]T q t )h m [l] (q u q 2 , q 3 , q 4 ) . (31) 

i=i 

where 7 4m [l](gi, q 2 , q 3 , q±) is the standard 4-mass scalar box integral. Previous arguments, 
based on the structure of the generalized form factor with two CPO insertions, imply that 
this should be the complete result. One may nevertheless check the relevant 2-particle cuts 
are correctly reproduced. Thus, this is the complete momentum-space form b 2 of b 2 . As 
in the case of the disconnected contribution, it is not difficult to Fourier-transform this 
expression and reeover the standard position-space form of b 2 : 

^ 2 = ~2 2 2 T ■ (32) 

x 12 x 23" ( '34 x 41 



26 Some UV regularization is assumed. 

27 To this end one treats independently the momentum of each propagator and introduces a 5— function 
for each three-point vertex. Using an integral representation of these four 5— functions makes all integrals 
trivial. 
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This expression is annihilated by the special conformal generators; consequently, their mo- 
mentum space expression (jSJ) should annihilate the 4- mass box integral (|3ip . 

4.1.2 Next-to- leading order 

As in the case or regular scattering amplitudes, we may either proceed systematically through 
the maximal cut method^! or, in simple cases, bypass some steps analyze directly a spanning 
set of color-dressed cuts. In the case at hand - the correlation function of four CPO-s - 
we may use the fact that we expect that it should be UV-finitd^l Thus, all contributing 
Feynman integrals should have at least five propagators. Ignoring disconnected and non-lPI 
contributions o, this implies that the contributions to the momentum space form of the 
connected part of the coefficient functions Oj, bj should be determined by the cuts in fig. 0] 
as well as the permutations of their external leg labels. As explained in § [2], each blob is 
either a color-dressed form factor or a color-dressed amplitude; depending on the specifics 
of operators, color-dressed form-factors of single-trace operators may be - up to the color 
factor - the same as color-ordered form factors. An example of such an operator is the CPO. 




Figure 4: Color-dressed cuts for connected part of the next-to-leading order 4-point correla- 
tor. There are 2 cuts of type (a), 4 cuts of type (b), 4 cuts of type (c) and 2 cuts of type 
(d). Other members of each set are related by relabeling. Each blob with an external line 
is a (generalized) color-dressed form factor; blobs with no external lines are color-dressed 
on-shell scattering amplitudes. 



The color algebra is nontrivial only for the cuts in figs. ID^a) andH(6). Writing the 4-point 
amplitude as 

A(h, l 2 , h, k) = Mi f*°* d f»°** + M 2 /«i«3 d/ a4»2d + ^ fai a idf a 2 a 3 d (33) 

28 Dcpending on the complexity of the correlator one may also use integral reduction strategy of [57] 
generalized to all-massive external legs, which reduces integrals to a basis at the level of maximal cuts. 

29 A divergence would signal the presence of an anomalous dimension for at least one of the operators, 
which contradicts the fact that the operators are chiral primaries. 

30 The cut in fig-SJa) also captures the non-lPI correction to the correlation function. The additional cut 
that is needed to determine the disconnected part may be interpreted as determining the next-to-leading 
order correction to the two-point function of two of the four operators. 
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and using the fact that the 2-point form factors are just Kronecker delta-functions in color 
space (with indices in the adjoint representation - e.g. the color factors of the two form 
factors at the left of the cut in fig. [4(a) are just S aia4 ) it follows that the cuts in figs. E(a) 
and H(o) are proportional to 

(M 2 - MO/^V^^Aa^ • (34) 

Further using the relation between Mi and color-ordered amplitudes it follows that the cuts 
in figs. 0(a) and 0(6) are 

%a) = ^(QiJiJ2)^(q2,kj2)^(q3,h,k)^(qA,k,k)M l ij3,kj2) (35) 

= F{qiAM)F{q2M,h)F{qzMM)F(q^M)M l uhM,h) (36) 

where A(/ l5 Z 3 , Z 4 , Z 2 ) is the color-ordered amplitude^! with the relevant external legsfEl This 
expression may also be justified using the photon decoupling identity. 

As in the leading order case, choosing as representatives of the four CPO-s the operators 
Ox = Tr[</> 12 12 ], 04 = Tr[0 34 34 ], 2 = Tr[0 13 13 ] and C 3 = Tr[</> 24 24 ] determines the 
contribution to the coefficient a 3 . It is not difficult to see that the cuts in figs. 0(o), (c) and 
(d) vanish identically because SO (6) invariance forbids propagators between form factors. 
The cut in fig. 0(a) is 



which implies that a 3 is given by the bow tie integral (see fig. 15(a)): 

r d A P dH i 

° 3 ~ J (2vr) 8 p 2 (p + qi ) 2 (p + qi + q A ) 2 l 2 (l + q 2 ) 2 (l + q 2 + g 3 ) 2 " 
This may be Fourier-transformed to position space with the result 



«3 



X A 1 X 



41^23 



d 4 x 



2 2 2 2 

x 01 x 02 x 03 x 04 

Permuting the labels of the operators we may similarly obtain a\ and a 2 . 

To determine the NLO correction to b 2 we choose as representatives of the four CPO-s 
the operators to be O x = Tr[</> 12 </> 23 ], 2 = Tr[0 14 24 ] O s = Tr[0 13 13 ] and £> 4 = Tr[0 24 34 ]. 
All cuts are nontrivial: 

(hh)(kk) (gi + g 4 ) 2 



(37) 



BTie(4,l|2,3) ; (38) 



(39) 



Ma) 



{hh)(hh) 
1 



(h + h) 



2Zi • h 



(hh)[hk] (kk)[l2k] 2 

Kkhnhu] (hh)[hh] 



m 
to^ 



U) 2 



(WW 

(hh){hh) (h 

(hh)[hk] 
(hh)(hl3)[hh][hh] ' 



(40) 
(41) 



31 This becomes the corresponding super-amplitude if one uses super-form factors. 

32 The coefficients Mj in the four-point amplitude ([55]) may be chosen to obey color/kinematics duality 
[35] on the unique off-shell leg of this amplitude. Thus, the corresponding cuts in figs. [4(a) and[4(&) obey it 
as well. 
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Figure 5: 2-loop integrals: BTie(4, 1|2, 3), DB(4, 1|2,3), TriP(4, 1|2, 3) and TriB(4|l|2, 3). 



Using the methods developed for the construction of scattering amplitudes of fundamental 
fields (e.g. matching onto an ansatz, or simply expressing the cuts in terms of momenta and 
inspecting the result) it is not difficult to find (the integral representation of) a function that 
has these cuts: 



£ (4) E ft) [tei + <?2) 2 DB(1, 2|3, 4) + ( qi + g 4 ) 2 DB(4, 1|2, 3) 



i=l 

+ gi 2 TriP(l|2, 3, 4) + g 2 TriP(2|3, 4, 1) + g 2 TriP(3|4, 1, 2) + g*IHP(4|l, 2, 3) 

- TriB(l|2]3, 4) - TriB(2|l|3, 4) - TriB(2|3|4, 1) - TriB(3|2|4, 1) 

- TriB(3|4|l, 2) - TriB(4|3|l, 2) - TriB(4|l|2, 3) - TriB(l|4|2, 3)1 (42) 

Here DB(a,b\c,d) is the double-box integral with legs a and b on one one-loop box sub- 
integral, TriP(a|&, c, d) is the triangle-pentagon integral with external leg a on the triangle 
point and TriB(a|6|c, d) is the triangle-box integral with external leg a at the 4-point vertex 
and external leg b on the triangle point; these integrals are shown in figs. [5J5), (c) and (d), 
respectively. 

The expression for 62 in eq. ( l4"2]) can be Fourier-transformed to position space; we im- 
mediately recover the expression for the integral representation of the coefficient 62 found in 
[51]. The easiest way to carry out the Fourier-transform is to (1) treat the momentum of each 
propagator as independent and introduce the appropriate 5— functions enforcing momentum 
conservation at each vertex (2) use an integral parametrization of the 5— functions to carry 
out the momentum integrals and (3) at this stage the Fourier-transform becomes trivial and 
identifies the integration parameters for the S— functions of vertices with external lines with 
the position of the operators. The resulting integrals (over the position of the vertices with 
no external lines) may be further reduced using a strategy outlined in the Appendix A of 
and we find the expected result [52]; we collect some of the details in Appendix O 



4.1.3 Brief comments on supersymmetric methods 

While the calculation in the previous subsection was carried out in components, it is not 
difficult to streamline it by making use of the super-form factor of the chiral stress tensor 
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multiplet. In such an approach it is not necessary to choose representatives of the CPOs; 
rather, products of the harmonic variables appearing in each form factor play the role of the 
various 6— functions appearing in eq. (j27|) . In the following we will continue constructing 
component correlation functions; in this subsection however we will illustrate the super- 
symmetric methods and recover the LO and NLO correction to the four-point correlation 
function of chiral stress tensor multiplets. 

The quadruple cut in fig. 3(b) is given bvF^l 

(T( qi )T(q 2 )T(q 3 )T(q4))\^ = ft) f[ ^ ^ / d%d%d%d% (43) 

5^( Tl+ - hrj[l] +;1 + Wl^^MlU - ^[1]_ ;2 ) 
£ (4) (72+ - kvM+fi + Z37/[2] +;3 )5 (4) (^[2]_ ;2 - l 3 r][2}_. 3 ) 
5 (4) (7s+ - W3] +;3 + ^[3] +; 4)5 (4) (^[3]_ ;3 - kv&Ud 
5 (4) (7 4+ - ^[4] +;4 + /ir/[4] +; i)5 (4) (^[4]_ ;4 - kvlM-.i) 

Carrying out the 77 integrals (either directly or though the methods described in pTI] ) 
we find the cut of a 4-mass box integral multiplied by the product of harmonic variables 
(12)(23)(34)(41) and 5 (4) (7i+)5 (4) (7 2+ )5 (4) (7 3 +)5 (4) (74 + ), implying that only the CPO com- 
ponents have non- vanishing four-point function [5]. 

Through similar manipulations and using the higher-point form factors ffTTj) it is not 
difficult to compute the cuts in fig. Ht up to a factor of £( 4 )(7 + i)5( 4 )(7 +2 )#( 4 ) (7+3)5^ (7+4), 
cuts (a) and (d) are: 

Cfc B) = (14) (23) + (12) (23) (34) (41) ( f + ^ + (13) (32) (24) (41) ( f + (44) 



= (12)(23)(34)(41) 
4 



(/1 + /3) 2 V /V /V /V / (/i + / 2 ) 2 

-(<?!+ <72) 2 + (/2-/5) 2 + (/2 + / 3 ) 2 , (/I 



+ 



(/ 2 + / 3 ) 2 (/2-/ 5 ) 2 (/2 + / 3 ) 2 (/2-/ 4 ) 2 

q\ -(gi + g 4 ) 2 + (/i + / 2 ) 2 + (/2-/4) 2 



(/i + / 2 ) 2 (/ 2 -/ 5 ) 2 (/i + / 2 ) 2 (/ 2 -/ 4 ) 2 J ' (45) 

We recognize the various momentum- dependent factors as the component cuts found previ- 
ously. Together with the cuts obtained by permuting the labels of the external lines these 
expressions can be used to recover the results of the § 14.1.21 



4.1.4 Higher-point correlation function of BPS operators 

The calculation of the correlation function of four scalar BPS operators may be easily ex- 
tended to the correlator of any number of operators. The main difference is that, unlike 

33 We use the standard convention p — > — p i-> \p) — > —\p), \p] —> \p]. The Grassmann coordinates rj[i] and 
77 [7] have different harmonic coordinates since they correspond to different multiplets - T]\x\°L j = iias^^a 
- where x labels the operator to which r\ belongs to. 
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the four-point correlator, the R-symmetry constraints are more difficult to write in compact 
form. Similarly to eq. ( 127)) . one is to identify the 5*0(6) singlets in the product of n 20'- 
dimensional representations of SO (6). By picking suitable combinations of operators and 
evaluating their correlation functions we may then extract the coefficient of each individual 
50(6) singlet. 




+ k + 1 i + k 



% + k + 1 



i + 1 



i - 1 



T [*] 



T M 



Figure 6: General two-loop integrals: I f , Iij , I2' ■ The integral Iif includes the numerator 

M 



factor s| = (Ej=o<?i+;) 5 



The cuts that need to be evaluated are similar to the cuts in fig. H] except that each form 
factor is replaced with a generalized form factor with arbitrary number of operators subject 
to the constraint that the sum of the number of operators is fixed to n. It is not difficult 
to see that the building blocks are identified by eqs. fl37j) . (JlO]) . fill) and they lead to the 
integrals shown in fig. El 

Based on the calculation of the 4-point correlation function one can see that integrals 
of the type Io[ enter only in the partly-connected (in R-symmetry space) components of 
the correlator and integrals of the type I]j and I2I appear in the connected components. 
Analyzing the cuts we find that the analog of the 6j coefficients in eq. (1271) . i.e. the coefficient 
function for the maximally connected 50(6) singlet, defined by the tensor structure 

n 

n<w> ( 4e ) 

8=1 

is given by 

n n—2 n n—2 

i=l k=0 i=l k=l 

where the integral Ii[ has a numerator factor sf^ = (q\ ) 2 with g) fc ' = X^y=o9*+j anc ^ 
first index on b denotes the fact that the index contraction is analogous to that of 62- 

While kinematic restrictions forbid a parity-odd part (i.e. containing Levi-Civita tensors) 
in the four-point correlation function, one may wonder whether such terms should be present 
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in eq. (]47[) . An analysis of cuts similar to those in figs. H^a), WLp) andH](c) shows that the 
integral topologies detected by these cuts do not exhibit a parity-odd component. Cuts 
involving two three-point form factors have two components, in which either one of the form 
factors is MHV while the other is MHV. While each component has a nontrivial parity-odd 
part, it cancels in the sum. Thus, to this order, no parity-odd part should be present in 
eq. flED- 

While it is not difficult (and perhaps useful for specific applications) to Fourier transform 
this expression to position space using the details in Appendix [C] and check its conformal 
invariance, the result is rather complicated. It is more instructive to carry out the Fourier 
transform under the assumption that we also take the null limit of the resulting expression 
and make contact with the n-point 1-loop MHV amplitude thus reproducing the results of 
[5J. Indeed, it is not difficult to see that the integrals Iq[ and are not proportional 
to Iir=i x ii+i an d thus are not sufficiently singular to contribute to this limit. The Fourier 
transform of integrals of the type Iij is (cf. eqs. (IC.16P and f)C.14p ) 

;n^V^)iif ] (47) 

3 

X l+k,i-l X 'i,i+k+l ~ X 'i,i+k X 'i-l,i+k+l + ®( X l,l+l) f 



I rf 2 / ,yi 2 ,y»2 ,-y» 2 ,-y»2 

Llj=l J A i,O x i+fc,0 i+fc+l,O x i-l,0 

for k 1 and k ^ n — 1 and subleading for these two cases. Assembling &i imax , the Fourier- 
transform of fci imax and after a trivial change of summation indices we recover, as expected, 
eq (4.11) of [5j. It is important to note that all contact term integrals in momentum space 
- i.e. integrals with cancelled propagators - do not contribute in the null separation limit 
of the Fourier-transform of the momentum space correlator. This might have been expected 
since at this order any contact term lies along a path between two null-separated operators 
and such terms have been argued [3] to be subleading. 



4.2 Correlators of BPS and non-BPS operators 

The methods described in previous sections can equally well be applied to correlation func- 
tions involving non-BPS operators, such as a general twist-2 operator. As discussed in § 12.21 
since these operators are nontrivial sums of monomials with coupling constant-dependent 
coefficients, the strategy for the evaluation of their correlation function is to evaluate the 
correlator of generic monomials and subsequently sum these components with the appropriate 
coefficients. We will illustrate here this construction with the calculation of the three-point 
function of one non-BPS twist-2 operator and two BPS scalar operators. We will then gen- 
eralize this calculation to an arbitrary number of scalar BPS operators as well as to two 
twist-2 operators. Compared to the correlation function of BPS operators, the additional 
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form factors which are needed are special cases of the general form factor found in § 13.41 



4.2.1 3-point correlators with one twist-2 operator 

For operators with general R-symmetry indices, 

Ox = Ti(D x + <f) h Dl~ x (f) Jl ) 2 = Ti(<f) h <f) j2 ) e> 3 = Tr(0 /3 Ja ) (48) 

the R-index structure of the correlator is 

(0,0,0s) = c? S[*6}?6* } + $s[*6%% } . (49) 

Depending on the choice of operators, both, one or none of the structures above are allowed. 
Whenever both are allowed the coefficient functions C\ and c 2 are related by the transforma- 
tion x -H- S — x. Here we will choose a particular non-BPS representative in which the two 
scalar fields are the same, I\ = J\\ in this case the factorization occurs on a term by term 
basis and its momentum/position dependence is given by (cf + cf). 

/ 13 




- 92 

Figure 7: Triple cut determining the leading order three-point function in eq. ( l49l) . The 
black dot denotes the non-BPS operator while the gray dots denote BPS operators. 



It is convenient to define 

S(a, b, x) = (a-f- x (b-) x + (b~f- x (a-) x , (50) 

where the upper index "—"denotes a projection onto a null direction (e.g. a~ = a + a 3 , 
etc). The leading order correlation function is determined by the triple cut in fig. [7] to be 

cf (0) + # 0) = / fl T^fe - h - k)S\q 2 + h- h)5\ gi + h + kf-^f (51) 

The Fourier transform to position space yields the expected result which may then be used to 
assemble the correlation function of the full twist-2 spin-S 1 operator with two BPS operators 
by summing them with the coefficients in eq. (125]) . 
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Figure 8: Cut topologies which determine the next-to-leading order correction to the three- 
point function in eq. (I49p . As in fig. the black dot denotes the non-BPS operator while 
the gray dots denote BPS operators. 



To construct the next-to-leading order correction to this correlation function we may 
proceed with the maximal cut method with the expectation that we will need to consider 
maximal and next-to-maximal cuts. At this order however the calculations are sufficiently 
simple to allow us to proceed directly to next-to-maximal and next-to-next-to-maximal cuts; 
up to the interchange of the two BPS operators (or, alternatively, up to the transformation 
Qi "H- q<2) they are shown in fig. |HJ These cuts are not difficult to evaluate using the form 
factors described in § [31 They are: 

c -»> = s{hM Wm- s(hM w^- (53) 

Due to the presence of two three-field form factors, the last two cuts receive contributions 
from the two configurations in which either one of them is MHV and the other is MHV. We 
find: 



c 



fig. Etc) 



+ 




S(h, k + h, x) - S(h + k, h, xj) 

\S(h + k,h,x) - S(h,l 2 + h,x)j + S(h + h,h,x) 
q\ 1 1 ql 



+ s{h > k + /3 ' x) ( (k + 1 A + w ~ Jhhf + wrw 1 1 (54) 
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%■ Erf) 



S(l x ,l 2 ,x) 



4 



+ 



+ 



4 



(k + U) 2 (k + h) 2 (h + h) 2 (l2 + h) 2 (h + U) 2 (h + h) 2 



+ 



4 



4 



(h + h) 2 (h-h) 2 (l2 + k) 2 (h-k) 2 (l2 + h) 2 (h-h) 2 
1111 



(k + k) 2 (k + k) 2 (k + k) 2 (k-hY 



(55) 



Inspecting these expressions and accounting for the symmetries of the cuts it is not difficult 
to find the (integrand of the) momentum space correlation function: 



+ c; 



■(1) 



S(h,l 2 ,x 




, 2 91 V\ 




+ % >\ 


+ ql 
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2 91 



+ 



-S(k + / 2 , / 3 , x) + l -pS(h, k + Z 3 , x) 

In 



k + 



k +i 2 

In 



S(k,k + k,x) - S(h + k,k,x) 



(S(h + k,k,x) - S(k,l 2 + k,x 




(56) 



As usual, each diagram stands for the product of scalar propagators corresponding to the 
graph, the product of the momentum conserving delta function at each vertex and an integral 
over all internal momenta with the numerator factor explicitly shown. 

While the spin of the operators can be arbitrarily large and thus the numerator factors can 
have arbitrarily high powers of loop momenta (cf. eq. ( |50|) ). Lorentz invariance guarantees 
that of all the integrals present in eq. ( 15 6 p only those with an explicit one-loop bubble 
sub-integral are divergent in the UV; this divergence is related to the fact that the non- 
BPS operator we consider acquires a nontrivial anomalous dimension (and also mixes with 
other operators). Thus, these integrals require regularization; we may simply promote to 
.D-dimensions the result of the four-dimensional calculation. The test of the validity of this 
procedure is the emergence of the correct anomalous dimension contribution to the correlator. 

It is not difficult to evaluate the one-loop bubble subintegral in dimensional regularization; 
from the Zj~ and (gi + /i)~ dependence of the result one can reconstruct the mixing matrix 
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of the monomials 0\ = O^g'^ 3 in ( )48l) . Diagonalizing it one finds the correct one-loop 
anomalous dimension - proportional to the harmonic number h(S) - and the corresponding 
eigenvectors, as well as the anomalous dimensions and eigenvectors of the descendants of 
twist-2 operators of lower spin. After taking the appropriate linear combinations ( 125]) 
of cf + c^ 1 ' to isolate an operator with definite anomalous dimension, multiplication of 
the full (LO + NLO) correlator by the appropriate renormalization factor (TTU]) will render it 
finite as the regulator is removed. To complete the calculation of the three-point function 
it is necessary to also include the tree-level correlation function fl5Tl) summed against the 
two- loop mixing coefficients [89] . 

It is relatively easy to check that the expression fl56|) vanishes in the appropriate limits. 
In particular, as S, x — > the non-BPS operator becomes BPS and the three-point function 
should vanish identically [92]. Indeed, in this limit the coefficients of the last two integrals 
vanish identically and thus the correlator is finite; after using momentum conservation the 
remainder is proportional to 




(57) 



We may similarly construct the correlation function of the spin-5 1 descendants of the BPS 
operator by noticing that it is given by the sum 



(58 



where Cf are the binomial coefficients; summing cf^ + cf W over integer values of x between 
and S the divergent integrals cancel out and the remainder is again proportional to ([5 



4.2.2 Correlators of one non-BPS twist-2 operator and {n — 1) BPS 

The calculation of the correlator of one twist-2 and two BPS operators can be generalized 
to correlators with any number of BPS operators. As in the case of correlators of only BPS 
operators, the cut calculation is almost identical, but in this case we need to evaluate the 
two additional cuts shown in fig. [9j These cut are very similar to the cuts in figs. M^c) and 
E^d), respectively, in which the intermediate state is restricted to only gluons and scalars. 
Similarly to the discussion in § I4.1.4[ the R-index structure is more complicated than that 

34 Thus, at this order the regulator dependence is completely accounted for by promoting the integration 
measure to be _D-dimensional. 

35 Upon Fourier-transform to position space this multiplication also restores the combination 71 In ^^^J,^ 23 ^ 
which is required by conformal invariance. 
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in eq. (H9"j) and it is given by all the SO (6) singlets that appear in the (non-symmetrized) 
tensor product of the representations of the BPS and non-BPS operators. We will focus 
on the completely connected part of the correlator and, as in the case of the three-point 
function, we will choose a non-BPS twist-2 operator with two identical scalar fields. 



</:> v , %-l Qi+1 , , %-l 




In ' x qk+l ' " q k+1 



Figure 9: Additional cuts necessary to determine the correlator of one twist-2 scalar oper- 
ator (with momentum qi) and an arbitrary number of dimension-2 scalar BPS operators. 
Momentum labels are slightly changed compared to those in fig. [HJ 



When non-vanishing, the leading-order connected component of the momentum space 
correlator of one spin-S twist-2 scalar operator and (n — 1) scalar BPS operators of dimension 
2 is 



C^ = J\{d%5{ qi + k-l l+l )^^ 



(59) 



where the external leg q\ corresponds to the non-BPS operator. Denoting by l\ and l 2 the 
momenta of the internal lines immediately adjacent to the external line qi - corresponding 
to the source-field of the non-BPS operator - in the integrals in which this line is not shown 
explicitly in fig. [TOj the NLO correction to this correlation function is: 



n n—2 n n—2 

[k] 
L 2i 



(7« = S{h,l 2 ,x) £I>i* ] -£I> 

\i=l k=l i=2 k=l 

^S(h + t 2 , t 3 , x) + ^S(h, t 2 + t 3 , x) 

3 (s{pi,P2 +P3,x) ~ S(p 1 +p 2 ,p 3 ,x)) hi 

(S(p 1 +p 2 ,P3,x)-S{p 1 ,p2+p 3 ,x)) h-,i ■ (60) 
P2 v J 

The integrals appearing in this expression and their internal momenta are defined in fig. [TUJ 
A potentially non-vanishing odd part present in various contributions to the cuts in fig. [9] 



P2 

^ Pi + P2 
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i + k i + k + 1 i + k i + k + l i + l i + l 




(a) (b) (c) (d) 



Figure 10: Integrals Iij , I2I , ^3,1, entering the n-point correlator of one non-BPS scalar 
twist-2 operator and any number of dimension-2 BPS operators. The first two integrals are 
the same as Iif and I 2 f in fig. [6j 



cancels in the complete expression of the cut. This is consistent with the absence of a 
parity-odd component of all the other cuts. 

With this expression in hand we may test whether the non-BPS nature of one operator 
affects the relation between the null limit of the correlator and the MHV n-point amplitude; 
general arguments in a specific (non-standard) regularization scheme [3] suggest that the 
relation is unaffected. The similarity of the contact term integral Fjf with the integral 
with the same name entering the correlation function of BPS operators (1471) and the fact 
that the Fourier-transform of the latter is subleading in the null implies that the integral 
l 2l here will also yield only subleading terms; as in that case, it is not difficult to check 
this explicitly. Inspecting ( 160]) and comparing it with ( 14"T|) we notice that, up to a factor 
S(h,l2,x) = (li) x (^2~) x + ^i) x (^2) S ~ x the integrals are identical. This factor is present in 
the leading order correlator (1591) and this cancels out in normalization. Of the other integrals 
in (1601) . I3 and I4 could also give contributions of the same order. If however the null limit 
is taken in the presence of a finite dimensional regulator, these integrals are subleading by 

and x^, respectively, and thus drop out as well [5j[3]@ We therefore explicitly confirm 
that the null limit of the correlator of one non-BPS and (n — 1) BPS twist-2 operators is the 
same as that of the correlator of n BPS twist-2 operators and the same as the MHV gluon 
amplitude. We also note that, as in the case of correlators of BPS operators, momentum 
space contact term integrals do not contribute to this limit. 

4.2.3 The correlation function of two twist-2 and any number of BPS primaries 

Generalizing slightly the calculation it is not difficult to construct the correlator of two 
non-BPS twist-2 operators and any number of BPS operators. For general choices of BPS 
scalar operators and twist-2 operators one would have several possible choices of R-index 

36 If the null limit is taken after the correlator is expanded in small e at finite xf i+1 , the contribution of 
these integrals should factorize into a scheme-dependent coefficient function [5]. 
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structures; the coefficient of each one of them may be obtained either by choosing suitable 
representatives or by using supersymmetric methods. We will focus here on the structures 
that have a nonzero position-space tree-level term if the twist-2 operators are constructed 
from Tr[DS- x <p AB D x _(j) AB } and its conjug ateH There are two classes of R-index structures 
that can appear: those in which the two twist-2 operators are adjacent (i.e. their R-indices 
are contracted) and those in which the two twist-2 operators are not (i.e. their R-indices 
are not contracted but rather are contracted with the R-indices of BPS operators). 




(o) (b) 



Figure 11: Additional cuts needed to find the correlation function if two non-BPS operators 
and any number of BPS operators at NLO. Cut (a) is absent if the scalar fields in the two 
twist-2 operators are not conjugates of each other. 



Apart from the cuts already used to construct the correlators in the previous section, we 
need the information provided by the cuts in fig. fTTT a) in order to find the correlators in the 
first class. Without loss of generality we will assume that the two non-BPS operators are 
at locations 1 and 2; whenever there are only two internal lines attached to the non-BPS 
operators, we denote their momenta by l±, I2 and I2, h, respectively. The result for the NLO 
contribution to this correlator is: 

(n n— 2 n n—2 \ 

EE 1 *? 1 -EE 1 *? 1 
i=l k=l i=3 fc=l / 

37 I.e. the scalars in the two BPS operators are conjugates of each other. 
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n-2 , _ _ x 

+ E ( + t2 ' *3> x ) + Sr^i. *2 + *3, z) ) 5(Z 2 , h, y)I 2 f 1 

t— 9 \ 2 2 / 

S'^i + * 2 , * 3 , y) + f^fa, t 2 + t 3 , y)) z 2 , x)i 22 fel 
(<S'(pi,P2 + P3,y) - S(p 1 +p 2 ,p 3 ,y)^S(p 3 ,p 3 - q u x) I 3;2 



n-2 

+Ei 

k=2 


^2 


n— 2 

, V i 
+ 2J 

k=l 


f Pi 


p7 + 


Th , 
— 1 


P2 




Th + 


— i 


* 2 




! pr + 








! P2 + 


ft. 


p 2 





(5(pi,p 2 +p 3 ,x) - S(pi +J?2,J?3,a?))S , (Pl +P2,Pl +P2 - ?2,2/) Ia;i 

5(pi +p 2 ,P3,x) - S(p 1 ,p 2 +P3,x)^jS(p 1 ,p 1 - q 2 ,y) h-i . (61) 



Note that the position of the bubble integral in the line connecting the two non-BPS operators 
is ambiguous; we have chosen it such that the coefficients of and I 3;1 are the ones that 
appear in the anomalous dimension of the twist- 2 operators. Such an organization exists for 
any operator. 

It is not difficult to extend this result to an arbitrary number of BPS and non-BPS 
twist- 2 operators in a split configuration (i.e. with all adjacent BPS operators). The only 
difference compared to eq. ( 161]) is the proliferation of integrals I3 and I4: each twist-2 non- 
BPS operator contributes two such integrals in a pattern easily identifiable by comparing 
the past two terms in eq. ( IBTJj) and the last four lines of eq. (IBTj) . 

For the correlators in the second class, in which the non-BPS operators are not adjacent, 
we need to also evaluate the cut fig. fTTT b). in which both such operators enter via a three- 
point form factor. It turns out however that this cut does not reveal any new terms. Let us 
assume that the two non-BPS operators are at positions 1 and j and that, whenever there 
are only two internal lines attached to each of them, their momenta are li, Z 2 and lj, 
respectively. The NLO contribution to the correlator correlator is given by: 

(n n— 2 n n— 2 

EE 1 *?' 1 - E E 1 ^ 1 
1=1 k=l i=2,ij^j k=l 
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+Ei 
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P 2 
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P2 
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p 2 





pr5 , (t 1 ,t 2 + t 3 ,y))5'(Z 1 ,Z 2 , : r)I 2 f ] 

J 2 / 



(<S'(Pl,P2+P3,^) - S'CPl +P2,P3,a:))5'(^,^+l,y) l3;l 

^5(pi + p 2 ,p 3 ,x) - S(pi,p 2 +p 3 ,x)jS(lj,l j+1 ,y) I 4; i 
S(pi,p 2 + p 3 ,y) - S(pi +P2,P3,y])S(h,l 2 ,x) h-j 
S(pi +p 2 ,p-3,y) - S(pi,p 2 +p 3 ,y)\S(li,l 2 ,x) U-j . (62) 



In both cases the potential parity-odd components cancel (up to total derivatives) in the 
complete expression of cuts. As in the case of correlators with a single non-BPS operator, 
it is not difficult to Fourier-transform eqs. ( 161]) and (|62|) in the limit in which the operator 
insertions are null-separated. As in the case of a single non-BPS operator, the null limit is 
to be taken at finite dimensional regulator. 



5 On effective actions and energy flow correlators 

As discussed in the introduction, an interesting application of the techniques we described is 
the construction of effective actions in background fields — in particular in non-dynamical 
gravitational background — as well as the construction of correlation functions describing 
the energy and charge flow in scattering processes of gauge-singlet states [6j [7] . Since the 
relevant correlators are somewhat different in the two cases we will discuss them separately. 



5.1 Effective actions 



The linearized coupling of a field theory with a background gravitational field is universal: 

S = So + J d d x bTT^ + ... (63) 

where Sq is the flat space action, is the departure of the background metric from that of 
the flat Minkowski space, is the stress tensor and the ellipsis stand for terms non-linear in 
h. The effective action for the gravitational field is obtained by integrating out all the fields 
except for W . We may interpret this as the evaluation of the scattering amplitudes of the 
background gravitons h^ u . As discussed in § P2J up to terms in which two gravitons emerge 
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from the same vertex, this is the same as the evaluation of the stress tensor correlation 
functions; that is, up to contact terms that may perhaps be constructed on the basis of 
general coordinate invariance and other symmetries, the effective action for the background 
gravitational field is nothing but the collection of the stress tensor correlation functions. To 
construct them in our approach we need to understand how to extract the stress tensor form 
factors. We will describe this in the next section. We will however leave for the future the 
complete evaluation of effective actions, including the determination of contact terms either 
by requiring that all Ward identities are satisfied or by finding and using the multi-stress 
tensor form factors. 

5.2 The MHV stress tensor form factor 

Thus, both for the evaluation of effective actions as well as for the evaluation of energy 
and charge correlators it is necessary to evaluate correlation functions involving the stress 
tensor perhaps together with other operators, perhaps with additional constraints in the 
structure of the contributing diagrams. They may - in principle - be constructed by a 
careful application of supersymmetry transformations on correlation functions in which the 
stress tensors are replaced by dimension-2 BPS operators and separation of the contributions 
of the descendants of lower levels in the supersymmetry multiplet. Extraction of the stress 
tensor component is, however, not obvious since it is necessary to separate from the coefficient 
of 9 2 9 2 all the contribution of the supersymmetry descendants of the operators appearing on 
lower levels in the multiplet. 

To evaluate such correlators in our approach it is necessary to know the form factors 
and generalized form factors of the stress tensor; their MHV components are part of the 
MHV super-form factor of the complete stress tensor multiplet constructed in [SI] , which is 
a natural generalization of the super-form factor of the chiral stress tensor multiplet: 

11 1 

q = ny-T^^-E^^ ( 64 ) 

i=l ' ' ' i i i 

n 1 

= n/-^ 5(4) ^-E A ^ (8) ^-E A ^)- ( 6s ) 

i=l ' ' ' i i 

This super-form factor and its generalizations may be used to construct super-correlation 
functions from which the desired components can then be extracted. Alternatively, one 
construct directly the component correlation function by first extracting the relevant form 
factors. 

A naive extraction of the component form factor from fl65|) as the coefficient of 7 + cr M 7-_7 + <7 iy 7_ 
leads to expressions which violate the expected conformal Ward identities related to the trace- 
lessness and conservation of the stress tensor. This is a reflection of the difficulty mentioned 
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above, that apart from the stress tensor, the coefficient of 9 2 9 2 also contains descendants of 
the operators at lower levels in the supersymmetry multiplet. It turns out that the stress 
tensor form factor {i.e. the coefficient of 9 2 9 2 from which the descendant contribution is 
separated) may be extracted as 

(0|TH1, ...,n) = ^jd 8 1 e ABCD lA ia» lBlc v lD G . (66) 

Conservation and tracelessness follow from the fact that e ABCD ^AlB = 0, which is a conse- 
quence of the Grassmann nature of 7 and of the antisymmetry of the Lorentz contraction. 
We have checked that this expression agrees with a Feynman diagram and BCFW-based 
evaluation of this form factor. 

As discussed in section |2j to construct correlation functions though generalized unitarity 
we also need generalized form factors with more operator insertions. In the case of the stress 
tensor they are particularly important to consider because the derivatives present in the 
stress tensor lead to contact terms which are missed if only regular form factors are used; it 
is possible that such contact terms - which are important for the construction of effective 
actions - can be determined by requiring that correlators exhibit conformal invariance. 



5.3 Energy and charge flow correlators 

Energy and charge correlation functions [HI IH2J EH ESj are particular examples of event 
shapes, which are important observables in QCD. For some some initial state X and observ- 
able e(N) that depends on the final state, they are defined as weighted cross sections: 

1 dcr t 



'tot 



J2da x ^ N 5(e-e(N)) (67) 



0tot de 

and describe the distribution of e(N) in the final state of some scattering process while being 
insensitive to the particular composition of this state. Popular choices for e is the energy or 
(some) charge flowing in some direction, e.g. 

N 

£(n)\ki, . . . , /cat) = ^(cos^ — cos 6^)5(0 — fa) k® . . . , fcjy) (68) 
i=i 

and similarly for some charge Q. The angles (9i,4>i) define the unit vector = ki/\ki\ with 
respect to some reference direction (9, fa). The operator £ admits presentations [961 E] i n 
terms of the stress tensor: 



d 3 n5(n 2 - l)w(n)£(n) 

/+00 
dtn'T^ifaRf?) . (69) 
-00 
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Using the expression of the differential cross section in terms of squared S-matrix elements 
one casts [Ml I9~7] eq. ( 157|) into an expression whose building blocks are special correlation 
functions of the operators £. In the definition above £ is understood to contain a projection 
onto states at t = +00. 

In a conformal field theory the initial state X is created by the action on the in vacuum 
state of some local gauge invariant operator: 

\X) = 0|O) in (70) 

Then, the analog of the total weighted cross section is 

1 da tot in (0\O^(e- Jd 3 nw(n)£(n))O\0) in 



(Ttot de (0|CW|0) in (71) 

This expression may be refined by assigning a different 5— function factor to each S-matrix 
element Using an integral representation of the 5— function it is easy to see that the expec- 
tation value (ITT]) is defined by the more general correlators [9] 



G(m, - - - n n , O) = (£(m) . . . £{n n )) = in(O | t | O)m " (72) 

The fact that the vacuum states are both in-states makes this correlator different than the 
standard ones which were interpreted as scattering amplitudes of source fields. Rather, 
as mentioned in [95], both the numerator and denominator have the interpretation of cut 
diagrams 0; we emphasize that these are regular unitarity cuts that have an imaginary part 
interpretation: in this case it is the imaginary part of the two-point function (0|O' O\0). The 
fact that £ contains a projection onto out-states implies that, in these unitarity cuts, the 
operators £ are to be inserted only on the cut lines; this is a reflection of the definition f )67|) 
which assigns the energy weight factors only on the final state \N). For the same reason, 
as well as because £ acts in the final state on widely separated particles, all operators £ 
commute with each other and thus no ordering needs to be specified. Singularities may 
nevertheless arise if some = fij and the corresponding operators act on the same particle. 
We will avoid such situations by keeping 7^ fij for all i 7^ j. 

5.4 On the evaluation of energy correlators at weak and strong coupling 

In contrast to correlation functions of stress tensors, energy correlators [HI [7] provide a class 
of observables which are insensitive to contact terms. Indeed, as discussed in detail in 5 15.31 



38 0ne may similarly define charge correlators, in which (some of) the energy detectors £ are replaced by 
the R-charge detectors 

/+00 
dtn 1 J l {t,rn l ) . 
-00 

39 This is the same as the cut diagram interpretation of inclusive cross sections. 
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the stress tensors featuring in eq. 072p are connected only to cut lines and therefore, all 
potential contact terms are set to zero. Moreover, since the energy measurement is made at 
infinity it affects only one particle at a time and thus only the two-point form factor of the 
stress tensor is necessary. 

As described earlier in this section, the energy correlators in the state created by the 
operator O may be interpreted as the unitarity cut of the two-point function of O and 
dressed with energy- and orient at ion- dependent factors. The discussion in § [2] implies that 
this construction may be given an interpretation in terms of the form factors of O. Indeed, 
inserting complete sets of out-states we have 

in (O|0t0|O) in = ^ in (0|Ot|A0(Ar|e>|0) in , (73) 
x 

in {U\0^£{n x ) . . . £(n n )O\0) in = in^liV') • ..S(n n )\N)(N\O\0) in , (74) 

N,N' 

where the sums run over all physical Fock space states \X) of some fixed momentum; we 
recognize each factor as the form factor of O and its conjugate. 

To evaluate the right-hand side of ( 174]) one may use the original definition (I5B|) of S, 
superficially perhaps at the expense of not using supersymmetric methods to evaluate the 
sum over intermediate states. We may alternatively use the stress tensor super-form factor 
and carefully take the limit R — > oo limit at the end while enforcing the relevant on-shell 
conditions. As we shall see shortly, the limit R — > oo requires that the momentum of the 
stress tensor operator vanishes; in this case it turns out that the form factor is just the 
traceless part of p tJ, p u , where p is the momentum of one of the external fields and thus 
factorizes in the sum over intermediate states. One may show along the lines of [97] that, 
for x = (t, Rn), 

Urn H 2 j dt j ^EL^ELsM( p ^S< + >(pI) (rfn,p') e «" e -»"F( Pl ,p 2 , . . . ) 

The argument of the integral on the left-hand side of the equation above is the typical term 
in the momentum space correlation function corresponding to (1721) after the state sum was 
performed and before the limit R — > oo is taken. F represents the contribution of the two 
form factors and of the other stress tensors except for one; the contribution of the last stress 
tensor is just {Piriip^). Evaluating the angular integrals in the saddle-point approximation 
(justified by the large value of R) leads to the orientation of the space-like parts of p\ and 
P2 be the same as the external vector n. The integral in the time direction sets to zero the 
energy carried by stress tensor operator thus setting p\ = p\ and thus also sets equal the 
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norms of the space-like components of the momenta p\ and p 2 and thus also \X) = \X'). 
This identity therefore allows us to use supersymmetric methods to evaluate ([71]). 

In [7J it was suggested that the energy correlation function may be extracted by suitable 
analytic continuation of the time-ordered correlation function. It is in principle possible to 
use the form factor interpretation described above to evaluate them at strong coupling; the 
construction the form factors [98] in terms of a Y — system would be an essential ingredient. A 
difficulty arises however in the treatment and cancellation of infrared divergences: while the 
energy correlators are infrared-safe observables, form factors exhibit infrared divergences. In 
a fixed order calculation they cancel against real emission contributions - lower-loop higher- 
points contributions to form factors. It is not obvious how the cancellation proceeds at 
strong coupling; we leave this important problem for the future. 



5.5 Examples of energy correlators 

We will discuss examples of energy correlators for two-field operators. The simplest example 
involves a single operator £. At leading order the calculation is trivial; we need the explicit 
expression of the two-point form factors of the desired operator and we need to carry out 
the phase space integral. For definiteness let us consider the scalar BPS operator, also 
discussed in [7J. In this case the form factors are unity; assuming that the state created by 
the operator has 4- momentum (q,0), using the notation p = p/\p\ and using the symmetry 
under the interchange of the two cut lines as well as the on-shell conditions we find 

m (0\O^(n)O\0}^> = 2 J ^g- ^g- Ul 6W(n - pJPfc + fiWm + UJ 2 -q) 

2 f ufdu! . . 2 q 

-<J(2wi - q) = — ^ — (76) 



(2vr) 6 j 4ui ^ (2tt) 6 16 

d 3 pi d 3 p 2 



47T f ufdUi c/n 4:71 1 



(2nr J 4^1 6{2Wl ~ q) = j2^8 (7?) 
Taking the ratio (172"]) we find [7J 

<^)>bps = £ • (78) 

The spherical symmetry is, of course, due to the state created by C?bps = Tr[(0 12 ) 2 ] being 
spherically symmetric [7J. 

We may introduce an asymmetry by considering the state created by an e.g. twist-2 oper- 
ator. Since in this case the operator is a sum of terms, to reconstruct both the operator and 
its conjugate, we need to consider the mode general building block m(0\Ol s £(fii)O2 ! s,x\0)m 
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with the unspecified R-symmetry indices chosen such that the total R-charge vanishes. It 
evaluates to: 



(2tt) 3 2u;i (2it) 3 2u 



2 



o m»„ . ns r(l^)' " (79) 



(2?r) 6 16 \2J v ' \l + n 3 

(p s 1 + x p 2+ )(Pi+ v P v 2+ )^(pi + P2)<*(wi + w 2 - g) 



(27r) 3 2u;i (2tt) 3 2w 2 

1 1 /g\2S /■ 2 25 /l-cos0 



- - 



, (1 + cos 9Y b , (80) 

(2vr) 6 8 \2) J v 7 \l + cos9J ' V ; 

where we used that pi !+ = p^o +Pi,3- Summing over x and y the ratio of the two expressions 
above with coefficients given by (12^|) and using the orthogonality of Gegenbauer polynomials 
to evaluate the angular integral we find that the leading order of the energy one-point function 
in the state created by a non-BPS twist-2 operator is 

= ^ (S + \) CT (n 3 ) 2 • (81) 

— * 

For twist-2 operators expressed in terms of an arbitrary null vector i = (£ , €) the energy one- 
point function is obtained from the above by the replacement n 3 H- t ■ n. By appropriately 
changing the coefficients we may also construct the same one-point function in the states 
created by conformal descendants of twist-2 lower-spin operators. 

The correlator of two (or several) energy flow operators is an observable with nontrivial 
parameter dependence even in the case of a spherically-symmetric state. Let us briefly 
discuss the correlator {£ (Hi) £(H 2 )) bps- We find: 

in (0|Oj 3PS £:(n 1 )£:(ri 2 )OBPs|0) i ( n 0) 

(2^20,1 (2*)*L 2 " Pl)U26{2){fi2 ~ p2)63(fl + p2)6{Ul + U2 ~ q) 

TtAS ^2<^ (2) (^2 - p 2 )u 2 5 {2) (^ -P 2 )5 3 (pl + p 2 )5(uj l +u 2 - q) 
{2it) 6 2uji (2n) 6 2uj 2 

M d \Z (o d \Z ^ifl (2) ("2 - Pi)uJi5 {2) (n! -p 1 )5 3 (p 1 + p 2 )5(uj l + u 2 - q) 
{Zixy^Zui \lixy :> luj 2 

j^y 6 D (^\n x + n 2 ) + 8® (ft - n 2 )) ; (82) 
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together with (177)) it implies that 

(£(ni)5(ffe))g s = £ (5 (2) (^i + H 2 ) + 5 (2) (^i - fkj) • (83) 

Thus, this correlator vanishes identically as long as Hi 7^ ±H 2 - The same 5— function de- 
pendence on the fixed vectors Hi seems to be generic for the leading order two-£ correlator 
in the state created by any two-field operator. At higher orders this appears to change; 
below we will explore the next order - which is also the first nontrivial order - assuming 
that Hi 7^ ±n 2 . 




Figure 12: The four components of the cut contributing to the energy correlators at the 
NLO in the state created by a BPS two-field operator. The operators £(H) are placed on 
the cut propagators. 



The components of the cut contributing to energy correlators at the first nontrivial order 
(we will refer to it as the LO) in the state created by a BPS two-field operator are shown 
in fig. [T2J They may be evaluated directly using the one-loop two-point form factor and the 
tree-level three-point form factor of dimension-2 BPS operators or they may be obtained by 
constructing the cut of the momentum-space NLO two-point function. Using the discussion 
in sees. [2J [3] and H] it is easy to find that, up to a factor of the square of the momentum carried 
by the operator, the two-point function is given by the two-loop scalar bubble integraQ 
From the integral form of the two-point function is it easy to find the four components of 
the cut such that momentum labels are aligned. 

A calculation similar to the one leading to eq. ( 182)) shows that the contribution of the 
first two components is proportional to 6^ (Hi ± H2), which vanishes due to our assumption 
that Hi 7^ ±fi2- The last two components are the same up to relabeling of momenta; we 
will analyze only the first one and double it. The contribution of both energy operators are 
inserted on the same cut line is proportional to 8^ (Hi — H 2 ) and thus may also be ignored. 

40 This is consistent with the non-renormalization of the two-point function of BPS operators: the two-loop 
bubble integral is proportional to q~ 2 and therefore the momentum-space two-point function is a constant. 
Fourier-transforming to position space yields only contact terms, proportional to 5 4 (xi —X2). It is interesting 
to note that in this case the energy correlator is related to such contact terms. This is special to states created 
by BPS operators. 
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The remaining terms are: 



^OIOU^iK^^bpsIO)^ (84) 
2 f <Pp 1 d 3 p 2 d 3 p 3 (3) . 

1 1 

2(wiu; 2 - Pi ■ P2) 2(a; 2 w 3 - P2 • Pz) 
x ( Wl u; 2 (<y< 2 >(ni - Pi)5 (2) (^2 - P2) + 5 (2) (™2 - pO^fa - £,)) 
+ Wl u;3(5 {2) ^i -pi)5 (2) (n 2 - p 3 ) + 5 (2) (™2 -pi)^(ni - &)) 
+^ 3 (5 (2 H^i - P 2 )5 (2) (^2 - ]3 3 ) + 5 (2) (n 2 - p 2 )5 (2) (^i - Ps)) 

The six integrals are essentially the same; we will discuss only the ones on the first line 
of the parenthesis. The integral over the angular parts of p 1>2 sets pi i2 = toi j2 n a> b with 
a 7^ b G {1,2}. Then the p 3 integral sets p 3 = — {uJin a + W2^&) and w 3 = |p 3 |. Last, we 
can do the integral over \p 2 \ = u 2 ; the S— function fixes u 2 and introduces a Jacobian factor 
J — u 3 /(q — u>i(l — fix ■ n 2 )). Most factors cancel and we are left with (for < n x • n 2 < 1) 



2 1 r/ 2 WidW! 



32(27r) 9 1 - f?i • ra 2 J (q- - ni • ra 2 )) 2 
g 2 (1 - ni • ra 2 ) - (1 + f?i ■ r? 2 ) ln(2/(l + fh ■ n 2 )) 

16(2tt) 9 (l-fii-fi2) 3 (l + ni-fi 2 ) 

Similarly, the integrals on the second and third lines of the parenthesis are 



55) 



T Q 2 -2(1 - gi ■ n 2 ) + (3 + gi ■ n 2 ) ln(2/(l + fk ■ n 2 )) 

1,3 16(2vr) 9 (l-ni.^l + fJi.^) 1 j 

/ 2 ,3 = il,2 (87) 

Combining them, accounting for the component in fig. [T27 d) and normalizing^] by dividing 
by dZZD leads tc@ 

{t{ni)t{n 2 )) B p S - — — - ^ — — . (88) 

2 (27r) 4 (1 - m ■ n 2 ) z (l + m ■ n 2 ) 

It should not be difficult to proceed to higher orders or to evaluate the energy correlators 
in states created by other operators. For example, at the next order one may start from 
the presentation of the two-point function in terms of the three-loop ladder and crossed- 
ladder propagator integrals, construct their three- and four-particle cuts and dress them with 



41 Since the LO correlator in(0|Og PS 5(ni)£(n2)C?BPs|0)i^ the cut of the NLO 2-point function contributes 
only to NLO energy correlator. 

42 This logarithmic dependence is the same as in QCD while the coefficient of the logarithm is different. 
We thank Lance Dixon for emphasizing this point to us. 
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Figure 13: Plot of the two-energy correlator in the state created by a twist- 2 spin-S* operator 
for S = (BPS), S = 4 and S = 8 and for S = 16, S = 32 and S = 64 as a function of 
ni • r?2, for £ • ni = = £ • n 2 - The energy scale is set to q — 1 and may be restored to 
arbitrary values by conformal invariance. 



energy factors. In evaluating the remaining integrals one sound observe the cancellation of 
IR divergences between the contribution of these two cuts. 

We may also easily evaluate the two-energy correlator to first nontrivial order in the state 
created by a non-BPS twist-2 spin-S 1 operator. The relevant dressed cut is: 



i 

d 3 pi 



o (pi +P2+ P3)0{ui + uj 2 + w 3 - q) 



d 3 p 2 



(2tt) 3 2u;i (2tt) 3 2cc; 2 (2tt) 3 2u;3 
q 2 



16(wiw 2 -pi •p 2 )(w 2 W3 

Pl+ + P2+ 
8p 2 +(^1^2 -Pi -V2) 

+S{ Pl - 



P2 ■ Ps) 



-S(p 1 +p 2 ,P3,x)S(pi,p 2 +P3,y) 



S(pi+P2,P3,x)(S(p 1 ,p 2 +P3,y) - S(pi +P2,P3,V)) 



x ( oj^iS 2 ^ - pi)5 2 (n 2 - p 2 ) + & 2 {ni - p 2 )5 2 (n 2 - 
+ ujxU^ 2 ^ - pi)<5 2 (n 2 - p 3 ) + 5 2 (rai - p 3 )5 2 (n 2 - pi)) 
+ u-iU^S 2 ^ - p 2 )5 2 (n 2 - p 3 ) + 5 2 (f?i - p 3 )5 2 (n 2 - p 2 ))\ 



(89) 



The evaluation of integrals proceeds following the same steps as for (£(ni)£(ri 2 ))BPS- The 
resulting expressions, which are now functions of four variables - S, ni-n 2 , i-n\ and £-n 2 where 
i defines the null direction in the twist-2 spin-S" operators - are quite lengthy. Multiplication 
by 8/(2tt) 2 (2S + l)/g 25 yields the desired two -energy correlator, (£(^i)£(n 2 )) 2) s- 

Its expression simplifies somewhat by choosing the vector i defining the light-like direction 
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to be orthogonal to the two unit vectors Hi^, i-e for I • Hi = = £ ■ n 2 . In fig. [13] we show 
the graph of {£(Hi)£(n 2 ))2,s f° r 5 = 0,4, 8 and for 5 = 16,32,64 for this choice of £, as a 
function of Hi • H 2 - Similarly to the two-energy correlation in the state created by a BPS 
operator (recovered here as the 5 = case), it has a simple pole at Hi ■ H 2 = —1 and a 
pole of order (25 + 4) (5 ^ and even) at Hi ■ H 2 = +1. The numerator exhibits the same 
logarithmic dependence as in the case of a state created by BPS operators, i.e. it depends on 
ln(2/(l + fii -H?)). Unlike that case however, the coefficient of the logarithm is a polynomial 
of degree (25 + 2) in Hi ■ n 2 ; the numerator also contains a polynomial of degree (25 + 3) 
in n\ ■ n 2 which factorizes as (1 — n\ ■ n 2 )P 2 s+2{Hi ■ H 2 ) and thus lowers the strength of the 
pole at Hi ■ H 2 — 1). These polynomial factors are a consequence of the momentum factors 
present in the operators. Indeed, setting 5 = we recover the expression in eq. (|88|) . By 



comparing figs. 13(a) and 13(b) we also note that the details of the operator with definite 
anomalous dimension are important for the two-energy correlator in a large spin state and 
are not important at low spin. Indeed, it is easy to see that the oscillatory behavior of 
Gegenbauer polynomials is not reflected by (£(n 1 )£(n 2 )) 2jS & t low spin, fig. 13(a), while it 



is obvious at larger spins, fig. 13(b) 



The same general features persist also if £ • Hi 7^ 7^ £ ■ n 2 . These scalar products do not 
affect the pole structure of the logarithmic terms of the two-energy correlator; they appear 
only in the various polynomial factors due to the momentum dependence of the operator 
creating the state. 



6 Outlook 



Correlation functions of local gauge-invariant operators are an important class of observables 
in both conformal and non-conformal field theories. In this paper we extended the gener- 
alized unitarity method to the calculation of momentum space correlation functions; apart 
from the tree-level amplitudes of fundamental fields, correlation functions of gauge-invariant 
operators Oi are determined by their on-shell form factors and generalized (multi-operator) 
form factors. This approach makes use of the methods developed for efficient construction 
of on-shell scattering amplitudes and exposes properties of some of the contributions to the 
integrand of momentum space correlation functions. In particular, we have argued for the 
presence of color/kinematics duality for all internal propagators which do not begin or end 
at an operator insertion up to contact terms that collapse at least one line connected to an 
operator. Using this approach we have recovered the known expression [511 [52] for the four- 
point function of scalar BPS operators of dimension 2, we have found the n— point functions 
of such operators at the next-to-leading order and discussed the limit in which the operator 
insertion points are light-like separated. As expected [5], we found that in this limit the 
correlator reduces to the one-loop n— point MHV amplitudes of Af = 4 theory. The methods 
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described in this paper are not restricted to correlators of BPS operators. To illustrate this 
we have computed the next-to-leading order correction to correlation functions of one and 
two non-BPS twist- 2 spin-S* operators and any number of scalar BPS operators of dimension 
2 and described the features of the correlator of m twist- 2 and n BPS operators in a split 
configuration. We have moreover discussed the M = 4 realization of the energy correlators 
and evaluated the one-point and two-point correlators to leading nontrivial order in states 
created by the action of the chiral primary and twist-2 non-BPS operators. It is not difficult 
to extend these calculations to higher orders. It would be very interesting to understand the 
relation between the standard approach we followed and the analytic continuation prescrip- 
tion suggested in [7|. It would also be very interesting to understand how to carry out the 
same calculation at strong coupling, perhaps by making use of the available techniques [HE] 
for the evaluation of form factors. 

The use of generalized unitarity for the construction of correlation functions follows 
closely the construction of scattering amplitudes and thus it produces momentum space 
correlation functions in an algorithmic way. It would be interesting to formulate generalized 
unitarity directly in position space; unlike its momentum space version, such a formulation 
would likely manifestly preserve conformal symmetry. It is possible that it will bring further 
simplifications to the Lagrangian insertion method. 

As suggested and illustrated in [51 [801 EE] position space n-point super-correlation func- 
tions of the chiral stress tensor multiplet are, in the null-separation limit, proportional to 
n-point superamplitudes. We have reproduced this result to leading nontrivial order using 
the momentum space correlation functions constructed here; we have also explicitly shown 
that the presence of non-BPS operators does not spoil this conclusion. 

The calculation of general correlation functions at large 't Hooft coupling remains an 
interesting open problem. The use of form factors, for which Y-system-based methods have 
been developed [98J, together with some variant of unitarity is a possible approach, as dis- 
cussed in § 15.41 Worldsheet form factors, recently discussed in [HH] , should provide an alter- 
native approach to the same quantities. To this end it is necessary to find the worldsheet 
form factors of the vertex operators dual to gauge theory local operators and use them to 
construct correlation functions of vertex operators using worldsheet generalized unitarity. 

Our approach suggests that, at least in certain situations, recursive methods could be 
used to construct the higher-loop integrand of momentum space correlation functions. To 
this end a necessary ingredient is the existence of a BCFW shift of operator momenta for 
which the tree-level form factors and generalized form factors are well-defined at infinity. 
Most operators explicitly discussed in this paper have this property, as can be seen from 
sees. 13.21 I3.3[ 13.41 and Appendix IB. 31 Clearly, if such a recursion relation exists it will 
determine simultaneously scattering amplitudes, form factors of operators, as well as their 
correlation functions. While we believe that a recursion relation exists - if only because, 
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for correlation functions of BPS operators at generic positions, recursion relations have been 
argued to exist at strong coupling [2H E2] - the nonplanar nature of correlation functions 
appears to be a (hopefully temporary) obstruction to formulating such a recursion relation 
in weak coupling perturbation theory due to the non-uniqueness of the parametrization of 
non-planar integrands. 

As we illustrated briefly in § I4.1.3[ supersymmetric methods initially developed for the 
construction of form factors [M| 153] apply naturally to the calculation of correlation functions 
and are quite efficient if all operators are annihilated by the same supersymmetry generators. 
In this case the presence of Grassmann 5— functions in form factors and scattering ampli- 
tudes simplifies the evaluation of sums over internal states and reduce^ them to simply an 
exercise in solving linear algebraic equations [91] HUUj . While for non-BPS operators such 
convenient presentations of form factors do not exist, it is possible to construct a coherent 
state representation for them once a component is found by some means. Adopting such a 
representation |101j for both the form factors of BPS operators amplitudes seems likely to 
simplify the evaluation of state sums for correlation functions involving non-BPS operators. 
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Appendix A: Harmonic superspace conventions 

The harmonic superspace [82J provides a bridge between a non-chiral M = 4 off-shell su- 
perspace and the chiral on-shell momentum superspace used to specify external states of 
scattering amplitudes and form factors. By introducing harmonic variables (u^ a , u^ a ,) val- 
ued in the coset 



an SU(4) index is split into indices transforming under a local 577(2) x SU{2)' x U(l) 



SU(A) 



(A.l) 



577(2) x 577(2)' x U(l) 



subgroup: 



6 A ^ 9 + ' 





(A.2) 



This holds both in TV = 4 sYM theory as well as in theories with fewer supercharges. 
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By gauge-fixing the local symmetries one obtains a parametrization of the coset in which 
the SU(4:) symmetry is broken. Keeping the harmonic variables generic allows one to carry 
out the projections (1A.2I) while maintaining an unbroken SU (4) symmetry. 

As an element of SU(4) the 4x4 matrix {u Aa ,u Aa ,) is unitary and thus its 4 x 2 blocks 
have the following properties: 

ui a u+ b = 5 a b , u^u M , = 5^ u^'u+ b = ui a u- Ab , = 0, (A.3) 

where the bar stands for complex conjugation. Similarly, 

<*+ + = Si ■ (A.4) 

Unitarity also requires that the matrix have unit determinant: 

1 ABCD + ab + - c'df - _ -I ( \ K\ 

I 6 U Aa e U Bb U Cc' e U Dd> ~ 1 > {A.O) 

where we use the convention: 

<?* = <?■'* = l = -eu = -e vv . (A.6) 
Eq. (1A.5|) leads to the relation: 

2 e U Aa e U Bb --U_ € c > d ,U_ . (A.7J 

There exists one set of harmonic variables for each operator; rather than attaching another 
decoration to u, it is convenient to use the operator label (or the insertion point label) instead: 
u Aa ^ k~Aa- quantity which commonly appears in calculations manifestly preserving the 
SU(4) symmetry is 

_ l c ABCD h + c ab h + /+ c cdi+ _ lt.Cc' ZDd' /+ c cdi+ (\ o\ 

[LZ) - i e k Aa e k Bb L Cc e L Dd - -^k_ e c , d ,k_ l Cc e l Dd . (A.8) 

The on-shell superspace coordinates r] can be projected similarly to (1A.2j) : 

7] A ^ V a + = ufriA v- = u A _ a 'r, A . (A.9) 

and similarly project the M = 4 fields such that the on-shell superfield is invariant; this 
amounts to inserting a suitable number of identity matrices in the form (I A . 4 [) in the su- 
perfield in footnote Fourier-transforming the coordinates 77" leads to the projective 
(non-chiral) superspace of ref. |102] which may be obtained by dimensional reduction from 
a six-dimensional (1, 1) superspace. The superfield is then given by [102] : 

+v%r +2 + v 2 v a 'i;a> + \v 2 v 2 <f> 43 (A.10) 

where 77 = 77+, f\ is the Fourier conjugate of r?_, rf = \rfr\ a and fj 2 = \f\ a iff ' . 
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Appendix B: 



In this Appendix we summarize the calculation of the various form factors and generalized 
form factors used in our examples of calculation of correlation functions of local operators 
through the generalized unitarity method. The following shorthand will be convenient: 

Pij = ~(Pi+ Pi+i • • • + Pj) , Q~-j = % + Pij ■ (B.l) 

Furthermore, we will use e for spinor indices and e for the SU(2) indices of the harmonic 
variables and in addition to the definitions from Appendix [A] we introduce the following 
conventions: 

Afe^Aj = (ij) e\ = +1 . (B.2) 



B.l Tree-level (NMHV) generalized form factor of the stress two tensor mul- 
tiplets 

As may be easily seen using Feynman diagrams, the simplest component of the generalized 
form factor with two insertions of the chiral stress tensor multiplet involves two scalar oper- 
ators connected by a propagator. We may also interpret this as the result of an MHV vertex 
expansion; this particular generalized form factor arises from two MHV vertices, as shown 
in fig. [2(a); we will therefore refer to it as NMHV. 

In computing the NMHV form factor, one SU(4) index will be treated differently than the 
others. We will denote it by T, while we will sum over the indices A, B, and C. By breaking 
manifest SU (4) symmetry we can get a relatively compact expression while still having the 
possibility of making a convenient choice of T. For the derivation of the non-contact part of 
the color-dressed components of this form factor we will use the MHV vertex expansion, in 
which one sums up the contribution of MHV diagrams with the topologies shown in fig. [2(a). 
It is not difficult to see that 

-cNMHV 

1 - 



i=l j=i r =i \™e>%3 I \J -°6 ij I ij-i 



(B.4) 



where we have already combined the Grassmann 8— functions of the two MHV form factors 
to expose the overall super- momentum conservation constraint. The factor 8^ may also be 
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written as the four-fold antisymmetric product of its argument. Using two such factors we 
partial fraction the ratio of spinor products as 



-cNMHV 



1 



3!IT m=1 Km+l) 



-5 8 (l a + 7i + + 2+Tf + " E VmX m )S\qi +q~ 2 -pi---- p n ] 



m=l 



n i—2 



X 



i=l j=i 



X 



X 



\P\6 
A j A i-1 



+ 



+ 



+ 



A j+l A i 



(jQl) (i - lQi,) (QU) (jQj,) (QU + lQj,) (Qj,i) (QU + 1) 



(B.5) 

To simplify somewhat this expression we introduce a factor of (ij) / (ij) or similar except 
in the cases where this is factor takes the value 0/0 (perhaps due to restricted kinematic 
configuration). Using repeatedly Schouten's identity as well as the relation 



;b.6) 
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we find, after a lengthy calculation, that 



-pNMHV 



1 n 

3' FT (mm + l/ {l+a1 ^ + 2+a 72 > ~ 2l V m \ m )6\ gi + q 2 - Pl 

■ 1 lm=l \ ' "1" / m= i 



Pn) 



TABC 



n i~2 



E E (4 a (7i + Qi) - E^^yx 1 ^ - E^a k ) 

i=l ,7=1+1 r=i r=i 



AjAjAj (i-li) 



r=i 
n i— 3 



(B.7) 



r=i 

Af_ 1 Af_ 1 A^_ 1 (jj + 1) 



m 2 



+ 



\P \ 5 X A 

A j+l A j+l A j+l 



(i - li) 



E 



o b'c' } + 1 + 
oe L Ab' L Bc' 



2 

ii) 



The relation (iQjj) = (iQ}+ij) implies that we can shift the summation index and further 
simplify the expression of the form factor. In the end we find that the color-ordered MHV 
form factor is given by 



-cNMHV 



(B.8) 



3! rim=i( m > m + 1) 

2 n i-2 2 n 

EE E 4+EE^ 



Pn) 



m=l 



w=l i=l j=i+l 



w=l i=l 
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where we have used the definitions: 



A% =e TABC WQg) 2 - wUll + \Q7M +^VTr(r\Qf j \i])(wM + i) - J> Ar <rz)) 



r=i 

3 



(JJ + 1} 



(«S) a (^«) 2 0*Xv + i> 



Q/c' + + JTABC / 

*f = ~ 2(gff (~ W/c^ c <7^>(7^> - | e a ' d '<^^(7^74 + )(7i + 7i + ) 

+ ^w+ 6 e 6 c ( 7 ^ + i)5 a(i (7^ + 7| + ) - ^w+ 6 e b c ( 7 ^ + 2)5 ad ( 7 ^ + 7^ + ) 

Here the harmonic variable w takes two values corresponding the the harmonic variables of 
the two chiral stress-tensor multiplets. 

This analysis is sufficient for the construction of the generalized form factor of two CPO- 
s. For operators at higher levels in the chiral stress tensor multiplet it is possible that MHV 
generalized form factors exist. In this case it is necessary to supplement the terms above 
by the contribution of MHV diagrams in fig. ^ b). We leave the construction of the relevant 
MHV form factor and of F^ 4 ^ to the interested reader. 

B.2 MHV form factor of a scalar non-BPS operator 

The SU(2) sector of M = 4 sYM theory contains operators constructed out of two complex 
scalar fields whose two-point function vanished; it is perhaps the simplest sector which is 
closed under renormalization group flow. Operators in this sector are labeled by the number 
of fields of each type; operators with definite anomalous dimensions are linear combinations 
of single-trace monomials with fields ordered in different ways. In this Appendix we consider 
a generic single-trace monomial, C>A 1 B 1 ...A k B k = Trf^^ • ■ • 4>A k B k )', and construct its MHV@ 
super-form factor; we will find the result quoted in eq. (|23|) . By restricting the pairs (Ai, B,i) 
to take only two values one finds the building blocks of operators in the SU (2) sector. 

The MHV form factor can be computed easily when there are no gluons by just computing 
a single Feynman diagram; the minimum number of external legs is k. To add gluons, we 
will rely on a BCFW recursion relation. To preserve the MHV character of the form factor 
the added gluons must all have positive helicity. We will choose to shift the momenta of 
two scalar fields; because of this only three-point scalar-gluon amplitudes can appear in the 
recursion relation. 



44 While in its "standard" presentation all operators are constructed out of two holomorphic scalar fields, 
this sector may be embedded in the 50(6) R-symmetry in many different ways. 
45 Herc MHV stands for 2k Grassmann parameters - two for each scalar field. 
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Similarly to amplitudes, form factors carry color structure. It may be possible to gen- 
eralize the construction of [1U3] to color-dressed form factors (the generalization at least 
partly bypass the use of multi-peripheral color decomposition |104] since the presence of the 
operator makes the form factor potentially depend on other color tensors apart form the 
antisymmetric structure constants). Here we will find the color-ordered form factors from 
which the complete color-dressed one may be constructed. 

Let us focus on the form factor with external gluons and scalars and start by adding a 
gluon between two adjacent scalar legs, i and j, (in the sense that there is no other scalar 
between them); the rest of the external legs will be other scalars. We choose to shift the 
momenta of the scalars,: 

$) = \i)+z\j), (B.9) 
\j] = ■ (B.10) 

This ensures that there is no boundary term at infinity. Indeed, since the gluon has to be 
attached to either the leg i or the leg j, there will be one propagator dependent on z and 
the form factor goes to for z —> oo. The two BCFW contributions, Tj and Tj, are: 

rp = [n + m + lPn+l] 1 = + = (ij) 

[iP u+1 \ s ii+1 [a + i](i + ij)(u + i) (a + 1) (i + ij) ' 

rj= [W + i][< + i3]J_ = 0- (B12) 
[Pi+ijj] s i+ij 

The vanishing of the second contribution is a consequence of the fact that \j) = \j) and the 
existence of the necessary three-point amplitude requires that (i + 1, j) = i.e. the external 
legs cannot be in a generic kinematic configuration. 
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Proceeding to add one more gluon suggests that the expected expression for the color- 
ordered form factor with external scalars and gluons is 

F — C (B.13) 

m= a (kk+i) 

where % and j are the adjacent shifted legs, a labels the first scalar line before % (see figs. [T4l 
fT5|) and C denotes the contribution to the form factor that depends on the external lines 
j, j + 1, . . . a. To prove this conjecture we will show that the BCFW recursion relation 
preserves its form. 

Let us consider an arbitrary configuration of external lines such that there are at least two 
gluons between % and j and at least two gluons between a and % (other configurations must be 
considered separately). The possible contributions shown in figs. [TUfa) and ITffi b). in which 
there exists a tree-level amplitude with at least two gluons, vanish identically because the 
(desired) MHV nature of the form factor required all gluons to have positive helicities and 
thus requires that the amplitude factor vanishes identically (being related to am amplitude 
with a single negative helicity gluon). The BCFW contribution shown in fig. [Tffic) vanished 
because it requires that the momentum of the new gluon be related to the momentum of 
the leg j. The only non-zero diagrams are the ones shown in fig. [T51 Notice also that the 
computation does not depend on what exactly is on the right side of the graph between j 
and a; we will denote this part by C. The two contributions are: 

-pfig-HSU) _ {n + + lPg+i] 1 (a-P-n+i)(-h+ij) C 

[iiW] s u+ i (i - i-p ii+1 )(-p ii+x i + 2) ]T k -j a (kk + 1) nti 2 (" + !> 

C(ij) (i-li)(ai+l) 



UizKkk + i) (i-u + i) 

T &g.UMb) _ [A-iij -!][»- jj] 1 (a-Pi-ii) (-Pi-iij) C 

[A-ii?] n-ii (i - 2-A-iiX-A-iii + 1) iTkZikk + 1} ntii(" + 1) 

C{ij) (of-l><ti+l) 



lltJaikk+l) (i-U + 1) 

In these expressions we used (IB. 131) for the necessary form factors 
Adding the two contributions jf g ' ^^- a ^ and T^ g ' leads to 



m=a(kk + l) 



C (B.16) 



which is consistent with the expectation flB.13j) . It is not difficult to repeat the analysis 
above for the special cases when the number of gluons between legs a and % and/or between 
% and j is or 1 and confirm the expression (IB. 16j) for these cases as well. 
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The computations with fermion legs are very similar. Making use of the fact that 
( 1B.16P is not sensitive to the R-symmetry index structure of the operator O A x B x ,...,A k B k = 
Tt(0 j 4 1 b 1 • ■ ■ 4>A k B k ), the form factors with n external scalars, fermions and gluons may be 
combined with ( TB.i6j) into the MHV super-form factor 

Fo AxBr ,...,A h B k = 1 ^ (mm+V> ^ (nH 0jAi6iSi ) Sp (f[V a . b . ) > ( B - 17 ) 

llm=l\ > ^ / { oll6l ..., ah) 6 h } \i=l / \j=l / 

where 

H-aAbB = VAaVBb ~ VBaVAb + SabVBaVAb, (B.18) 

(S« l6l )° 7 = >ZX e fh (B.19) 

and the first sum runs over all over the sets {a±,bi ■ ■ ■ ,cik,bk} with particle labels valued 
between 1 and k and ordered as a\ < bi < a 2 < b 2 ■ ■ ■ 6fc_i < a k < 6 fc or 6 fc < a\ < b\ < a 2 < 
b 2 ■ ■ ■ bk-i < a k etc, as stated in § 13.31 

B.3 Tree-level MHV form factor of twist-2 operators 

Another set of operators that is closed under renormalization group flow forms the so-called 
SL(2) sector; these operators are linear combinations of monomials of the type 

0$£ D = Ti(D x + <p AB D s + - x (P CD ) , (B.20) 
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Figure 16: An example of a form factor with two external fermions and one external scalar 




Figure 17: The two classes of BCFW terms contributing to the form factor of twist- 2 oper- 
ators. 



as stated in eq. ( 1B.20j) . In this appendix we show that the tree- level form factor of these 
monomials are given by eq. fT26|) 



?4( _ sr^n \ c-1 a-1 / k \ x ( I \Sx 

T °tx° = Tr q ~<mV+l) E H — H — EE Eft Eft < B - 21 : 

llm=l\ ' / { a>bjC)(J } k=b l=d \r=l+l J \s=k+l / 

(b\g-p k \c) + (6|^ + io-|c) _ / (d|o-^|a) + (d|^+i(r~|a) _ ^ 

the first sum runs over all sets {a, b, c, d} where a < b < c < d or d < a < b < c etc. Fig 
summarizes the notation we use. 

We will focus on the case when all external lines are either scalars of gluons. The form 
factor with two external scalars and no gluons can be easily constructed by inspection: it is 
just the product of the " — " components of the momenta of the two scalars to the appropriate 
powers. 

To prove eqs. (|26|) . (IB.21j) we will use a variant of a BCFW shift. In choosing the shift we 
would like to avoid a boundary term arising from z — > oo region; this can be accomplished 
by the following shift that involves the momentum of a gluon and of the operator itself: 

\i) = \i) + z\cf>) : (B.22) 
^ = ^+ (B.23) 



57 



(a) 



(b) 



Figure 18: 



where \<p) is chosen to satisfy [z|er~|</>) = 0; such a choice is always possible. With this 
shift the BCFW recursion relation for the MHV form factor of 2 s'^ D has only the two 
contributions shown in fig. [TTJ: 

^.nzta) j^i-Pia] i (p-y , - )S - X _ M ( - v - r(v ~)s-* 

{a{ayilb) (p- a+P mp;r* (b.24) 



2p i (ai)(ib) 



(Pa) X (Pi+P^) S - X ■ (B.25) 



(a|^a | b) 



2Pi (ai}(ib) 
in deriving these expressions we used that 

(a0) {i\a~\i]{a(j)) {a\a 



(B.26) 



which follows from the Schouten identity and the choice of \<p) such that [i\a \<p) = 0. 
Adding the T' 1 fig -Q21>) anc j j^g-EH^) we g nc j a resu it j s indeed consistent with the general 
expression shown in eqs. (1261) and (IB. 211) : 

(^SfW^^) (B.27) 

" (ai) <»&> (ba) \ [Pa } {P * + Pb } 2pr + {Pa + p < ] [Pb ] 2p~ ) ' 

This expression was used in the evaluation of some of the cuts in § 14.21 It is also possible to 
check that this agrees with the Feynman diagram evaluation of (0 2 s' x D \4 > ab^ + ^cd) ■ It is 
also easy to see that in the limit S, x — > eq. (1B.27|) reduces, as it should, to the three-field 
form factor of the BPS scalar operator. 

Let us proceed to the general case, in which we add one gluon to a form factor which 
already has some arbitrary number of gluons present; the two non- vanishing contributions 
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to the color-ordered form factor in the ordering in which all gluons are adjacent are shown 
in fig. dHJ 

S—x 



g-EEElv) _N [i-Pia] i i 



E ^+E?r E?; 

k=i+l \ r=i+l / \s=k+l 



[a-Pia] P iaU b m =i+i{m,m+l){P ia i+l) 

(Pia\<r-pk\b) | (P ta p k+1 a-\b) ^ A | (prW^p-j C^j^±i^) 



6 \ S-x 



1 / - , -yx f -1 / (a|g-^|6) (q|fo+i<r-|&) 



for +*r £ p; + v V - ^ 



+ 1 (a|<r-|»](ti + i) v ITp~\ i T v 

Iti<m,m + 1> 2pr<«+l> fe ^A^ 7 

x / (a|Q fc |6) + (a|^ + ia-|6) _ { S (fi 2g) 



6-1 / fc \ 1 / ». \ -S 1 -^ 



?2 



2 ft 2 Pfc+i 
fig.dSfe) _ [^ + 1] 4 1 



[-P«+ii][?i + !][«'+ ^tf+i n';i +2 K m + !> (aPu+i)(Pu+ii + 2) 

S-x 



" ii+l 



s=i+2 



6 x S-x 



+ e U + a« + x>-V ( e A s ~* ( { ^g^ + - w 

fc=i+2 \ r=i+2 / \s=fc+l / Z ^ 



n„i<™,™+i> vt? V 2 k n„„,('«.'«+i) 2 p. r <«+i> 



Adding together the two contributions we find exactly the scalar-gluon component of eqs. 
and (IB. 211) . We described explicitly here the construction of the form factor with adjacent 
gluons; a very similar calculation yields the form factor with some gluons above and some 
gluons below the scalar line. Among the MHV form factors there also some with two external 
fermions; For the appropriate color ordering it combines with the scalar-gluon form factor 
and yield eqs. (|2SD and flR2Tj) . 
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Appendix C: On the Fourier-transform of the NLO four-point func- 
tion of CPO-s 



Let us begin by Fourier-transforming the double-box integral to position space. 

DB(1,2|3,4) = |j^^^(27r) 2 V(g 1 + / 1 + / 2 )5 4 (g 2 + / 3 -/ 2 )5 4 (g4 + /4 + /5) 

5 4 (<? 3 + k - h)S% -h- k)S\-h -h- h) (CI) 



3=1 
6 



gia;4-(g44-k+i5)gia;3-(g3+i6— '5)g«a;5-(i7— '3— i6)gi:i;6-(— ii— '3— '7) 



e iE i=1 -i-*j 16 /l2 / 25 j 35 J 34 J 46 J 56 (C.3) 

3=1 



where Jjj are standard position space propagator 

Fourier-transforming ( 1C.3I) to position space simply freezes the integrals over xi,X2,x$ and 
£4 to the location of the operator insertions. Momentum space numerator factors become 
differential operators acting on the integrand. 

Repeating the derivation above for the other integrals in eq. ( 142]) we find that, up to 
derivative and numerator factors, the integrals that can appear in their Fourier transform 
are 



^123 = J d 4 wl lw l 2w hw (C.5) 

X\234 = J d 4 wli w l2 w l3whw (C.6) 
#12;34 = J d 4 ud 4 vli u l2u Iuvhvhv , (C-7) 

The correlation functions depend only on {d\ + c?2) 2 i^4i ; 23 or, perhaps more symmetrically, 

((d l + d 2 ) 2 + (d 3 + d 4 ) 2 )H 41;23 ; (C.8) 



following the ideas of [90], in this appendix we spell out some details of the reduction of the 
relevant derivatives of the H integral. 



46 In Minkowski space the propagator has an additional factor of ' 
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This integral does not require any regularization as it is finite both in the UV and in 
the IR. Furthermore, the terms in which two derivatives act on the same coordinate reduce 
immediately to integrals of the type upon use of 

= 5\x l -x j ) ; (C.9) 

We therefore need to reduce only 

Following [90], we multiply and divide by {x\ — X/i) 2 {x 2 — X3) 2 and rewrite the numerator 
in terms of derivatives on the integrand: 

(xi - x 4 ) 2 (x 2 - x 3 ) 2 
(xj. - x 4 ) 2 (x 2 - x 3 ) 2 

= -2^-2- " d P4 ) 2 e^ + ^) ((d P2 - d P3 ) 2 e^ +i ^) (C.11) 

X 1A X 23 J i=1 \ Z7l J 

P1-P2+P3-P4 , n v 4 r(4)/ . . . \ /n 1 o\ 

We then integrate by parts (d Pl — d PA ) 2 setting aside the factor (p^Psipi +P4) 2 ) -1 011 which 
this combination of derivatives act trivially and find that the integrand becomes 

(R _ R )2 Pi -P2 +Ps -Pi 4(pi + Pa) 2 (pi -P2+P3- Pa) 2(P2 + P3) 2 ( 1 _1 

\O pi O p4 ) ^2^2 ^4^4 ^2^2 1 „2 ~"~ „2 



PlPl PiPl PlPl \Pl PI. 

(C.13) 



2 2 2 2 

P2 _ Pi , Ps _ P2 

2 2 1 2 2 ' 2 2 

PiPl \Pl PI Pi Pi 



The main difference between the first two and the last terms above is that in the former the 
denominator factor (pi + Pi) 2 cancels out. We will therefore treat these terms separately. 

Plugging the last term in the original integral (lC.12p . extracting the remaining derivatives 
as and adding to the result its image under the combined transformations X\ «-> x 2 and 
X3 -H- 24 it is not difficult to see that this term leads to Y^-type integrals multiplied by 
ratios of distances x 2 ,: 

(( X 3A - ^13)^134 - O24 - ^ 2 )^124 + (^34 - X Ia) Y 234 - (X 2 l3 - x\ 2 )Y lM ) . (C.14) 



x 14 x 23 



A useful relation is 



x 2 12 Y 123 = 4 / d A y [ TT ^ e *<*^)^<^)^<*>^)-^(\ + \) (C.15) 
J J 7_i ( 27r ) P1P2P3 \Pi P2J 



d 4 pid 4 p 2 e ipi - Xl ' i+ip2 - X23 f d 4 pi e iprXl3 f d 4 p 2 e ip:i ' X23 



(2tt)8 pfpj J (2vr)4 p\ J (2tt)* V \ 
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which may be proven by starting with a Fourier transform representation for Y123. 

Returning to the first two terms in (10. 131) we notice that upon restoring the factor of 
{pip\{Pi +P4) 2 ) -1 the mixed-momentum propagator cancels out. Acting with the remaining 
derivatives, (d P2 — d P3 ) 2 , we find 

(f) _ rl ,2f -4(Pl-P2+P3-Pl) 2_ (l_ 1 

V t/p 2 ^P3 ) I 4 2 2 4 

\ PlPWi P1P2PWA \Pl Pi 

—32 

(Pi ■ P2P3 -P4+Pi- Pm -Ps-Pi- P3P2 ■ Pa) ■ 



4 4 4 4 
PiPiPlPl 



Integrating this expression with the measure factor following from eq. ( 1C.12I) one can check 
that the result may be reorganized as 

222222/" 1 
X 13 X 24 ~ X 12 X 34 ~ X 14 X 23 / U d y]e ipi ' X1+ip2 ' X2+ip3 ' X3+ip4 ' X4 (C 16) 

2xj 4 xl 3 J p\p\pIp\ 

Combining eqs. (1C 16|) and (1C 14[) it is easy to see that 

— jpX • 82 + 83 • ^4)-^41, 23 _ ^1234 _ ^1234 _ ^1234 
-^14-^23 Il?>l24 I34I12 h3h4 

where 



+ £^41 — (^2,41 + £4,32 — £1,32 . (C.17) 



G 3;41 = ^ - ^ ■ (C.18) 
J31 -'34 

Using similar techniques, or as a limit of (1C.17j) by noticing that in the limit x\ — > £2 
the right-hand side of eq. (1C.16I) becomes identically zero, it is possible to show that 

d\H A1 , Vi = ^Y 1U . (C.19) 
1 34 

which was originally proven in |5lj by explicit integration. 
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